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The pressure distribution due to a concentrated load on 
a semi-infinite elastic body is given by the well-known 
Boussinesq solution for either the two-dimensional or 
three-dimensional case. We here investigate the effect on 
such a pressure distribution, taken at the depth h, of the 
presence at that depth of a slippery rigid bed (case (b)); 
of a perfectly rough rigid bed (case (c)); and of a perfectly 
flexible but inextensible thin layer embedded in the 


material (case (d)). Both pressure distributions in the two 
and three-dimensional problems are calculated for each 
case, Fig. 5 and Fig. 6. Several authors have already in- 
vestigated case (b) and case (c) in two dimensions (2), 
(3). Their results are in perfect accordance with ours. 
The author is indebted to Dr. A. Casagrande for suggesting 
this investigation as a contribution to the field of soil 
mechanics. 





GENERAL SOLUTION OF THE EQUATIONS OF 
ELASTICITY 


N the following paper, solutions of the equa- 
tions of elasticity are needed which give both 
the stresses and the displacements. Therefore it 
will be easier to start from the equations relating 
the displacements u, v, w. 
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v is the Poisson ratio of the material. 

It was established by Neuber' that a general 
solution of these equations may be found as 
follows. 

We call &o a scalar function of x, y, 2 satisfying 
the Laplace equation V?¢o=0. 

We call © a vector, such that each of its 


1 Neuber, “Ein neuer Ansatz ziir Liésung raumlicher 
Probleme der Elastizitatstheorie,” Zeits. f. angew. Math. 
und Mechanik 14, 4 (1934). 


Cartesian components ¢, ¢, @, satisfies the 
Laplace equation 


V°o.=0, V%,=0, V2¢,.=0. 


A very general solution of the equations of 
elasticity (1) can then be expressed as 


 teaesllii (0/dx) (go+xo.+9b,+26:) +4(1 — v) dz, 
— (0/dy) (dot xb. +¥by+2o:) +4(1 — v) dy, (2) 
= (0/02) (dot+xd. +¥oby+2¢:) +4(1 = voz} 


or in vectorial notation 


v 


u= —grad (¢@o+r-)+4(1—)®. 


The expression for © becomes in terms of the 
arbitrary functions above 


©=div u= —div grad ¢)—div grad r-@ 
+4(1—») div €. 


Since div grad ¢9=V2o0=0 

and div grad r-6=V*(xo,.+ yo,+2¢,) =2 div $, 

we find QO =2(1—2v) div $. (3) 
The vertical stress component o, and the 

horizontal shear 7,.7,2 are the only stress com- 
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ponents we shall need in the following theory. 
They are given by 


o, Ow v& Tre OW OU 


7G a 1-2» G ax az 





(4) 
Tyz OW Ov 


G ody as 
where G is the shear modulus of the material 
and »v the Poisson ratio. 

It is important to note that the factor (1 —2v) 


disappears in the value of o, when using the value 
of © given by (3), 


o,/2G=dw/dz+2v div &. (41) 


This shows that the value of the term vQ/(1—2v) 
when v= 4, 0=0, tends toward the limit 2» div ®. 


THE Two-DIMENSIONAL PROBLEM 


Let us consider first the two-dimensional 
problem where the load is concentrated on a 
line and has the value P per unit length (Fig. 1). 
We shall assume that the ground is elastic, of 
elasticity modulus E, and shall restrict ourselves 
to the case where the Poisson ratio is v=}, 
which means that the material of the ground is 
supposed to be incompressible. Four cases will 
be investigated. 


(a) The ground is infinitely extended and deep. The 
pressure distribution at the depth fA is given by the well- 
known Boussinesq solution. 

(b) The ground is infinitely extended but of finite depth 
h. It is resting on a rock base, perfectly rigid; no friction 
forces are supposed to act between the rock and the upper 
ground. They can slip with respect to each other without 
the slightest resistance (Fig. 2). This case has been calcu- 
lated by various authors.’ Their result, derived by a dif- 
ferent method, is found in proper accordance with ours. 

(c) The same case as in (b) but where perfect adhesion is 
supposed to exist between the upper ground layer and the 
rock surface. No slippage whatever is allowed to occur; the 
two materials are assumed to stick together perfectly. 
This case was investigated by Marguerre* for y,=0. Com- 
parison with his paper shows that there is no practical 
difference between the cases y=0 and »=}. The pressure 
distribution is calculated at the rock surface (depth h). 

(d) There is a horizontal infinite inextensible but per- 
fectly flexible layer at the depth h. The pressure distribu- 
tion on that layer due to P is calculated. This case corre- 
sponds to the problem of a clay substratum containing a 
thin horizontal sand layer. 
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In cases (a) and (b) the pressure distribution 
is found to be independent of the elasticity 
constants of the ground, while it depends on the 
Poisson ratio v for cases (d) and (c). Hence the re- 
sult for cases (a) and (b) is not affected by the 
restrictive assumption that v=}. 


Case (a) 


We take the z axis directed positive downward, 
the x axis being at the surface of the ground, and 
we assume that all the xz planes are identical. 
It is possible to find solutions of the elasticity 
equations for which all the variables are cosine 
or sine functions of x by the assumption: 


oo= (Ae**+ Be?) cos Ax, 


(5) 
o:=(Ce**+De-**) cos dx, 1$:=¢,=0, 


where A, B, C, D, X, are arbitrary constants. 
By substituting this into Eqs. (2), (4) and 
(4,), we obtain; 


Oz 07h 06. dg: 
— = —— 2 $-A(1—»)—, 
2G 02? 02? 02 








(6) 
Trz 20 0g0 0g: 


a= SR a Z 


G OxL dz 02 





-—(1 ~2»)0| 











PRESSURE DISTRIBUTION 


Since we assume incompressibility, v=} 


a,/2G= —0°oo/d2* —206,/02°+09./02, 
t22/G= —2(d/dx)[db0/d2+20¢,/d2]. 


At infinite depth we must have o,=7,,=0 and 
this is only possible if A =C=0. We then have 


tz:/G=2d sin \x[ — Bye? —zDde~* J. 


By introducing the condition that the surface 
carries only a normal load and no shear, 7:2, 
must be zero for z=0, hence B =0. The expression 
for o, then reduces to 


o,/2G= — Dye? (1+2A) cos Ax. 


This last relation shows that if the surface z=0 
is loaded with a normal pressure distribution, 


—o.0= G0 cos \x, the pressure at the depth 
z=h is 


—o,=Gq=qoe*"(1+Ad) cos dx. (8) 


Now an arbitrary surface loading po(x) may 
be represented by means of the Fourier integral 


as a superposition of sine loadings of various 
wave-lengths. 


(os) +co 
po(x) = (1/7) f dn f dtps(£) cos M(x—8). 


The pressure distribution g(x) at the depth h is 
derived from relation (8) and the above Fourier 
integral. 


1 co (os) 
g(x) =- f dn f dtpo()e-*(1-4hd) cos (x—E), 


us 


which can be written 
1 +00 eco 

q(x) =-{ pala | ddye—*"(1+hd) cos A(x— &). 
TT? co 0 


If the surface loading po(x) is concentrated 


within a small region (—e,+.«) such that the 
total load is 


P= f Pied 


The pressure distribution at a depth h due to 
that load is 
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q(x) = (P/x) [ e™(1+N) cos Axdh, 
; (9) 


q(x) = (P/xi) {et +a) cos a(x/h)da. 


As will be seen, this yields the well-known 
Boussinesq solution. 


Case (b) 


The soft ground is supposed to rest on a rigid 
rock base, no friction occurring at the surface of 
contact (Fig. 2). This case is identical with the 
symmetrical loading illustrated by Fig. 3. In 
order to proceed as in case (a) and use the 
Fourier integral, we first have to consider a 
sinusoidal loading which is symmetrical with 
respect to the plane at depth h. 

We take the x axis at the rock surface, the 
surface of the ground being at z=h and its 
symmetrical image at z= —h, (Fig. 3). A sym- 
metrical solution of the equations of elasticity is 
found by putting 


oo0=cosh Az cos Ax, ,=A sinh Az COs Ax. 


The application of formula (7) gives the value 
of the horizontal shear. 


Tzz/G=2x sin Ax[A sinh \4z+A2d cosh Xz j. 


Since no shear is acting at the boundaries z= +h, 
we have the condition 

sinh AA+hAA cosh Ah=0 
or A = -—sinh \h/h cosh hy. 
With this value of A the value of ¢, is given by 
Eq. (7). 
O; |: sinh Ad sinh Az sinh Ah cosh dz 


- —cosh rs| 
h cosh Ah dA cosh Ah 


XX? cos Ax. 


2G 
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This last relation shows that if a pressure 
go COS AX is acting at z= —h, the pressure g on 
the plane of symmetry z=0, which represents 
the rock surface, is 

2(sinh Ak+Ah cosh Ah) 


qq=——— —Go COS AX. 
sinh 2Ak+2\h 





1 


Proceeding as in case (a) by the use of the 
Fourier integral, we find the pressure distribution 
at z=0 due to a concentrated load P at the 
surface z= —h 


2P f” acosha+sinh a x 
p=— —_—§—— cos ada. 


th Jo sinh 2a+2a h 





(10) 


This checks with the formula found by Melan.? 
It gives the pressure p at the rock surface when 
no friction occurs as a function of the horizontal 


distance x from the vertical line of application of 
the load P. 


Case (c) 


This case considers the same material, soft 
ground lying on rigid rock, but here the upper 
soft material is supposed to stick perfectly to the 
rock surface. 

We take the x axis at the rock surface and the 
2 axis directed positive downward, as in case (b), 
(Fig. 2), and start from a solution which is 
sinusoidal along x 


oo= (Ae**+ Be-**) cos Ax, 

o:=(Ce*+ De-**) cos Xx, (11) 
o:=9,=0. 

We have three boundary conditions 


u=Oatz=0, w=Oatz=0, 7r,,=Oat z= —h. (12) 


We remember that we assume the soft material 
to be incompressible (v=4). With this value of 
the Poisson ratio, we get from formulae (2) 


and (7) 
u= sin \x[ Ae**+ Be—?+2Ce**+2De-** |, 
w= —cos Ax[_(Ae**— Be“) 
+ Ce**(Az — 1) —De-*#(Az +1) ], 


?L. N. G. Filon, Phil. Trans. Roy. Soc. London, 1903; 
E. Melan, Beton und Eisen 18 (1919). 
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t::/G=2d sin \x[ (Ae** — Be-**)d 
+2h(Ce* — De-*?)], 
a./2G=X cos d\x[_Ce**(1 —Az) 
— De~**(1+Az) —d(Ae**+ Be). 


(13) 


The three boundary conditions determine three 
of the constants A B C D in terms of the fourth 
one. Calculating the coefficients B C D in terms 
of A, we find that the value of o, at the rock 
surface (z=0) is 


a,/2G=2A2? cos \x(cosh Ah 
+hd sinh \h)/dAh cosh dh 
and at the surface of the ground (z= —h) it is 


o./2G=2A cos \x(cosh? Ah+ (hd)?)/Ah cosh Ah. 


This shows that if a pressure go cos Ax acts at 
the surface of the ground, a pressure 


g = (cosh Ah+Ah sinh \)/ (cosh? Ah+ (Ah)?) 


is acting at the surface of the rock. 

Using the Fourier integral as in the previous 
cases, we find the pressure distribution p at 
the surface of the rock due to a concentrated 
load P acting at the surface of the ground. 





p(x) =— cos a-da. 
th» cosh? a+a? h 


P f cosh a+asinha s 
0 


This formula coincides with Marguerre’s result 
for the case v=}. 


Case (d) 


The soft ground is here supposed to be 
infinitely deep and to contain at the depth / an 
inextensible but perfectly flexible thin layer to 
which the soft ground sticks perfectly so that 
only vertical motion at that depth is permissible. 

We take the x axis at depth / on the in- 
extensible layer and the z axis positive downward, 
and a sinusoidal solution in x as in case (c). The 
only difference with the preceding case is the 
boundary condition (12). The vertical displace- 
ment w at z=0 is not zero, but is related to the 
normal stress on that layer. To find this relation, 


3’ Marguerre, “Druckverteilung durch eine elastische 
Schicht auf starrer rauher Unterlage,” Ing. Archiv 2 
* (1931). 














PRESSURE DISTRIBUTION 


let us consider a sinusoidal solution applying to 
the infinitely deep material located below that 
layer. This is a case analogous to case (a). 
Since the material is infinitely deep, we must 
have A=C=0 and the horizontal displacement 
u deduced from formula (2) with »=3 is found 
to be 


u= sin \x[_Be~*+2De-*? ]. 


The condition that the horizontal displacement 
be zero on the inextensible layer (z=0) is 
B=0. 

The vertical displacement at z=0 is 


w= De~**(1+ 2) cos Ax 
and the normal stress 
o,/2G= — Dde~**(1+Az2) cos Ax. 


The relation between w and a, at the surface of 
the inextensible layer is 


o,/2G= —w. 


The boundary conditions for the upper soft 
ground are then 


(a) u=0 at z=0, 
(b) o,/2G=—hw, 2=0, 
(c) Tz2>= 0, z=—h. 


As in case (c) we may write the relation (13) and 
introduce these boundary conditions to find the 
value of B C D in terms of A. 

The value of the normal stress o, at the rock 
surface z=0 is found to be 


C: cosh Ah+Ah sinh Ah 
—=2)» “COS AX 
2 hye” 





and the value of o, at the ground surface is 


o, 2A 


2G Ar 





[hd(hd cosh hA—hd sinh hA—cosh hd) 
+cosh hA+hxd sinh hd]. 


From this we conclude that a pressure dis- 
tribution go cos Ax at the ground surface trans- 
mits a pressure on the inextensible layer 





ote hAChn/ (1 +h tanh hd) —1] 
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By using the Fourier integral as in the previous 
cases, we find the pressure distribution p trans- 
mitted to the inextensible layer when the surface 
carries a concentrated load P 





PrP e-* x 
p=— cos a-da. 
tho 1—al1—a/(1+etanh a) | h 


THE THREE-DIMENSIONAL PROBLEM 


The same cases (a) (b) (c) (d) are investigated 
if the load P is not concentrated on an infinite 
straight line but concentrated on a point, (Fig. 4). 
The stress distribution must then be axial- 
symmetrical around the load. Solutions of the 
problem may be found by using axial-sym- 
metrical potential functions with cylindrical 
coordinates 7, 2. 


oo= (Ae**+ Be?) Jo(Ar), 
o.= (Ce**+ De) Jo(Ar), 
oy= 0, :=0. 


We have a vertical displacement w and a radial 
displacement u. Eqs. (2) and (3) become 


u= —(0/dr)[¢o+2¢. ], 
w= —(0/dz)[bo+2o.]+4(1—v)¢.. 


If we call 7 the horizontal shear acting in the 
radial direction, we have equations entirely 
similar to Eqs. (4), (41). 


o,/2G=dw/dz+v@/(1—2y), 
1/G=0w/dr+du/dz. 


As before, we shall assume that v=}. These 
equations show that the problem of determining 





" a 





the constants A B C D is the same as in the two- 
dimensional case. For instance, if in the two- 
dimensional case to a load distribution at the 
surface g 9 cos Ax corresponded a pressure at 


depth h 
g=2(Ah)go cos Ax 


we may conclude that in the corresponding 
three-dimensional case to a load distribution 


goJo(Ar) at the surface corresponds a pressure at 
depth h 


g = g(Al)goJo(Ar) 


and we do not have to repeat the above calcula- 
tions to find the function g(\h). 


Case (a) 


The case of a concentrated load P acting on 
an infinitely deep ground is found from the dis- 
tribution goJo(Ar). At a depth h the vertical 
pressure p transmitted is found by using the 
solution of the two-dimensional problem g(h)) 


=e *(1+4+hd) 
p=e-"(1+d)goJo(dr). 


An arbitrary load g(r) could be represented by 
using the identity 


q(r) = [af Jo(Ap)Jo(Ar)q(p)dpdp. 
0 0 


This shows that the corresponding pressure dis- 
tribution at depth h is 


pir)= [e-M(1-+ni)JvOwynar [Ta0do)a(o)ode. 
0 0 


If the load is concentrated at the origin on a 
circle of radius « ) 


P=2rf q(p) pdp. 
0 


The corresponding pressure distribution at depth 
h is 


p(r) =(P, 2n) f he**(1 +A) To(Ar)dd 
0 


or 
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p(r)= (P/2eht) [ ae~*(1+a)Jo(ar/h)da. 


The same method may be used for all other 
cases (b), (c), (d). We find: 





Case (b) 
P © acosh a+sinh a r 
p(r)= { 2e—_—__—_—— —"Jo( a) daw 
2h? 5 sinh 2a+2a h 
Case (c) 








P © cosh a+asinh a r 
p(r)= f a Ju( a”) da. 
0 


2rh? cosh? a+a? h 


Case (d) 








P eo e 
pir)=— [a 
2rh? J, 1—al1—a/(i+atanha) ] 
rT 
x Ju a") day 
h 


NUMERICAL EVALUATION OF THE INFINITE 
INTEGRALS 


The problem reduces to the evaluation of 
either 


(P/rh) f g(a) cos (ax/h)da 
0 
for the two-dimensional problem, or 
(P/2nht) f ag(a)Jo(ar/h)da 
0 


in the three-dimensional problem. The function 
g(a) has the form 


g(a) =(1+a)e@ 
acosh a+sinh a 


sinh 2a+2a 


case (a), 





g(a)=2 case (b), 


cosh a+a sinh a 
g(a)= case (c), 
cosh? a+a? 





e-2 
a)= d). 
g(a) 1—al1—a/(1+atanh a) | senathans 





We start from the two identities 
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Fic. 5. Pressure distribution p/(2P/mh) in the two-dimensional problem. 








Case (a) Boussinesq distribution. Maximum value of pressure p =2P/rh. 
Case (b) — — — Slippery rigid bed. Maximum value of pressure p = 1.441 -2P/zh. 
Case (c) Rough rigid bed. Maximum value of pressure p = 1.291 -2P/ah. 


Case (d) —-—-— Inextensible flexible layer. Maximum value of pressure p =0.935 -2P/rh. 
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Fic. 6. Pressure distribution p/(3P/27h) in the three-dimensional problem. 





Case (a) Boussinesq distribution. Maximum value of pressure p =3P/2xh?. 
Case (b) — — — Slippery rigid bed. Maximum value of pressure p = 1.711 -3P/2xh?. 
Case (c) Rough rigid bed. Maximum value of pressure p = 1.557 -3P/2xh?. 





Case (d) —-—-— Inextensible flexible layer. Maximum value of pressure p =0.942 -3P/2xh?. 
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few cos bxdx =a/(a?+5?), 


0 
f eeedulbxax = (a?+5?)-}, 
0 


Taking the derivatives with respect to a, we get 


d” a 
fovee cos bxdx=(—1)" ( ’ 
0 da” \a?+b? 


d” 1 
[ xrer*Julbxydx= (— 1)” ( ). 
0 da” \ (a*?+6?) , 


This shows that the integrals would be evaluated 
readily if the functions g(a) where all expressed 
as a sum of terms of the type a"e~**, where n is 
an integer. 














Case (a) 


In case (a) there is no difficulty in expressing 
these integrals. For the two-dimensional problem 





Case (b) 
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we get 


p(x) = (P/xb) { (1 +a)e* cos (ax/h)da, 


p(x) = (2P/mh)(1+-2?/h?)~. 


For the three-dimensional problem 
p(r)= (P/2eht) [ a(1+a)e~*Jo(ar/h)da, 
0 


p(r) =(3P/2eh?)[1+(r/h)2]-*?. 


These are both the well-known Boussinesq 
solutions. 
The two-dimensional distribution 


p(x) /(2P/mh) = (1+x?/h?)-? 


is represented by curve (a) in Fig. 5 and the 
three-dimensional distribution p(r)/(3P/27h?) by 
curve (a) in Fig. 6. 


The function g(a) may be represented with an error smaller than 1 percent by 


(2a cosh a+sinh a)/ (sinh 2a+2a)=2(1+a)e~*— (1+2a)e-** — 5. late **, 


The pressure distribution for the two-dimensional case is 


2 2 0.5 


p(x) = 





rhL(1+x2/h?)? (14+(x/2heP 





1 —6(x/4h)?+(x/4h)! 
059 | 
L1+(x/4h)?}! 


The dimensionless quantity p(x)/(2P/zh) is plotted in Fig. 5 by curve (b). 


For the three-dimensional case we have 


3P 2 
Ha 


0.25 








= —0.0 
rh Cit+(r/h)? 2 [1+(r/2h)?}*? 





- —3(r/4h)?+ | 
7 [1+(r/4h)2]9?2 j 


The ratio p(r)/(3P/27h?) is represented in Fig. 6 by curve (b). 


Case (c) 


The function g(a) may be represented with an error smaller than 1 percent by 


(cosh a+a sinh a)/ (cosh? a+a’) =(2acosh a+sinh a)/(sinh 2a+2a) —2.80ate—** — 56a‘e-* 8, 


The pressure distribution in the two-dimensional case is 


2 2 0.5 
p(x) = 





1—6(x/4h)?+(x/4h)! 
59 


1—10(x«/3h)?+5(x/3h)'! 
38 





- —0.0 
mhL[i+x?/h? PP [1+ (x«/2h)?? 


[1+(x/4h)?}* 





[1+(x«/3h)*]}> 


1—10(x/8.8h)?+5(x/8.8h)! 
—0.012 | 





[1+(x/8.8h)? ]} 
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| The ratio p(x)/(2P/7h) is represented in Fig. 5 by curve (c). The pressure distribution in the three- 








dimensional case is 

3PT 2 0.25 1—3(r/4h)?+2(r/4h)! 
p(r)= - —0.039 

2ehL[1+(r/h)2 2 [1+(r/2h)* 2 [1+(r/4h)2]92 


—0.15 





J —5(x/3h)?+(15/8)(x =) 
[1+ (x/3h)?}"/? 

The ratio p(r)/(3P/27h*) is represented in Fig. 6, by curve (c). 

Case (d) 

The function g(a) may be represented with an error less than 1 percent as follows: 

e~*/{1—al1—a/(1+e tanh a) ]} = (1+a)e-*—1.77a%e—*, 


The pressure distribution in the two-dimensional problem is 


2P 1 1—6(x/3h)?+(x/ 3h)! 
p(x)=—| ——___0.065 | 
mht (1+x?/h*)? [1+(x/3h)*}} 


The ratio p(x)/(2P/7h) is represented in Fig. 5 by curve (d). 
The pressure distribution in the three-dimensional problem is 








“USS 
rh (1+r2, h?)>/2 L1+(r/3h)?]*? 





3P [ 1 1—3(r/3h)?+3(r/3h)! 
p(r)= 058 | 


The ratio p(r)/(3P/27h?) is represented in Fig. 6 by curve (d). 
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Magnetic Properties and Orientation of Ferromagnetic Particles 
II. Pyrrhotite! 


C. W. Davis,? U. S. Bureau of Mines, Pittsburgh, Pennsylvania 
(Received June 21, 1935) 


The effect of heat treatment and of particle orientation on the magnetic properties of granular 
pyrrhotite are discussed and illustrated by curves and data. Heat treatment which causes little 
or no chemical change alters pyrrhotite so that increased induction and remanence and de- 
creased coercive force result. Particles aligned in a field parallel to that used in their magnetic 
measurement showed greater induction and remanence but less coercive force than randomly 
oriented material, while cross-oriented particles showed the lowest induction and remanence 


but the greatest coercive force. 





INTRODUCTION 
CONSIDERATION of recent observations 


on the effect of orientation on the measured 
magnetic properties of natural and of heat 
treated pyrrhotite shows that although the 


action with respect to coercive force is different - 


from that of magnetites it is nevertheless sus- 
ceptible to the same explanation that has been 
given for the behavior of the latter mineral.’ 


EXPERIMENTAL 


A massive specimen of magnetic, auriferous 
pyrrhotite from Haywood County, North Caro- 
lina, was crushed to pass a 100-mesh screen and 
was purified by repeated alternating and uni- 
directional field magnetic separation. The puri- 
fied sample was magnetically measured and was 
then heat treated at 250°C for 40 minutes to 
increase its magnetic induction so as to form a 
product better suited to measurement in our 
apparatus. 

The magnetic measurements were conducted 
as described in the previous article except that 
very small steps were used in the determination 
of coercive force and, further, that the sample 
was demagnetized by the recently described 
method in which a switch surrounded by a 
condenser is broken in such a way that a maxi- 
mum magnetic field inside a solenoid is reduced 
to zero by a damped oscillatory discharge.‘ 


! Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

2 Metallurgical division, U. S. Bureau of Mines. 

°C. W. Davis, Physics 6, 96 (1935). 

*C. W. Davis, Physics 6, 184 (1935). 


Particle orientation was performed as before 
except that cross-oriented samples were pre- 
pared by agitating and rotating the sample 
capsule with its axis perpendicular to the direc- 
tion of the field after each fractional addition 
of the sample was made. All results have been 
calculated to a packing density of 2.36 grams per 
cubic centimeter. 


DISCUSSION OF RESULTS 


The results show that there are certain 
characteristic differences between the effect of 
orientation on the measured magnetic properties 
of magnetite and those of pyrrhotite. This is in 
line with the difference in magnetic structure 
that has been attributed to the minerals by 
Weiss and his co-workers.°® 

They studied magnetic properties of thin 
plates cut in various crystallographic directions 
from crystals of magnetite and pyrrhotite and 
concluded that a few magnetite crystals ap- 
proached magnetic isotropism at field strengths 
of 400 oersteds or more but that most specimens 
had no greater symmetry than that of the 
rhombic system. They distinguished between 
two types of pyrrhotites: The normal, exempli- 
fied by crystals from Morro-Velho, Brazil, with 
a magnetic plane in which magnetization takes 
place more easily than in a direction perpendicu- 
lar to it; and the more common abnormal type 
in which the magnetic properties are poorly 
defined. They reported further that normal 


5 P. Weiss, J. de physique 4, 469 (1905). P. Weiss and 
J. Kunz, J. de physique 4, 847 (1905). Max Ziegler, 
Magnetic Properties of Pyrrhotite Crystals, Diss. (Ziirich, 
1915). V. Quittner, Ann. d. Physik 30, 289 (1909). 
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TABLE I. Analyses of natural and heat treated pyrrhotite. 








ANALYSES, PERCENT MOLECULAR 





OT RATIO 
DETER- 
SAMPLE DESCRIPTION Fe Ss MINED Fe Ss 
Natural, 100 mesh 56.97 36.75 6.28 10 11.23 
Heat treated, 100 mesh 57.09 36.05 6.86 10 11.00 


Heat treated, 325 mesh 56.95 36.47 6.58 10 11.15 








pyrrhotite may require a field of 13,400 oersteds 
for magnetic saturation. 

The heat treatment used to augment the 
magnetic properties of the North Carolina 
pyrrhotite caused superficial oxidation, as evi- 
denced by a resultant blue coloration, but pro- 
duced little chemical alteration, as is shown by 
the analyses in Table I. The ratio of iron to 
sulphur corresponds approximately to FeiSi, 
which falls within the composition range usually 
allotted to pyrrhotite. 

Magnetic properties of minus 100-mesh natural 
and heat treated samples are given in Figs. 1 and 
2. Fig. 1 shows that heating has raised the 
intensity of magnetization at H=3340 from the 
138 and 167 gauss of the random and oriented 
samples of natural mineral to the 286 and 353 
gauss of the random and oriented samples of the 
heat treated material. Fig. 2 shows that a similar 
change in remanent magnetization has resulted 
from heating, the increase ranging from 79 to 
123 gauss: for the random and from 82 to 158 
gauss for the oriented samples. The coercive 
forces, however, have undergone a reversal; the 
highest coercive force, 210 oersteds, is exhibited 
by the natural random sample, followed in order 
by natural oriented with 200 oersteds, heat 
treated random with 170, and, lowest, heat 
treated oriented with 135 oersteds (see Fig. 2). 

Since the change in chemical composition has 
been so slight, it would seem that heating must 
have affected the physical constitution, as 
discussed by Michel and Chaudron,® so that 
more perfect alignment of magnetic domains is 
possible. The decrease of coercive force suggests 
possible crystal growth on heating, inasmuch as 
Gottschalk’ has shown that there is a direct 


6A. Michel and G. Chaudron, Comptes rendus 198, 
1913 (1934). 

7V. H. Gottschalk, “The Coercive Force of Magnetite 
Powders.” Report of Investigations 3268, U. S. Bureau of 
Mines, p. 83 (1935); Physics 6, 127 (1935). 
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linear relation between coercive force and 
specific surface. 

Other tests showed that minus 325-mesh heat 
treated pyrrhotite gave results well suited for 
measurement with our equipment and different 
only in degree—not in kind—from results ob- 
tained with natural and with heat treated 
pyrrhotite samples of other particle size. 

Fig. 3 shows that for heat treated minus 
325-mesh pyrrhotite the magnetization (47/7) up 
to field strengths of 3340 oersteds is highest 
(404 gauss) for particle orientation parallel to 
the field, lower (292 gauss) for random orienta- 
tion, and lowest (266 gauss) for samples cross- 
oriented during capsule rotation. This 28 percent 
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Fic. 1. Magnetization curves for 100-mesh natural and 
heat treated pyrrhotite, different orientations. 
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Fic. 2. Hysteresis curves, upper branch near origin, for 
100-mesh natural and heat treated pyrrhotite, different 
orientations. 
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Fic. 3. Magnetization curves for 325-mesh heat treated 
pyrrhotite, different orientations. 


maximum difference in magnetization between 
the two first cases given above as contrasted 
with a zero difference at the same field strength 
in the case of the magnetites reported in a 
previous paper is matched by the magnetites 
when a comparison of their magnetization is 
made at a similar condition of unsaturation, a 
difference of 27.5 percent in their magnetization 
being obtained between oriented and random 
samples at a field strength of 42 oersteds. 

As the results for pyrrhotite are given for 
measurements made at about one-fourth the 
reported saturation value, it is very probable 
that in fields of about 13,400 oersteds the effect 
of orientation on the magnetic induction of 
pyrrhotite might disappear. 

Although the curves shown in Fig. 4 cannot 
represent absolute hysteresis values since satura- 
tion was not attained, the remanence as meas- 
ured doubtless shows the same relationship 
between different orientations that would ob- 
tain if absolute values were available; at any 
rate, the remanent magnetization differences 
follow the same general trend as shown by true 
hysteresis loops of magnetite samples. The same 
sample of pyrrhotite showed measured remanent 
magnetization (47/) as follows: Oriented parallel 
to the field, 208 gauss; oriented random, 136 
gauss; and cross-oriented, 126 gauss. 

The coercive forces determined from the data 
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Fic. 4. Hysteresis curves, upper branch near origin, for 
325-mesh heat treated pyrrhotite, different orientations. 


used in the construction of the curves in Fig. 4 
are also subject to the unfortunate limitation of 
our measuring equipment with a maximum field 
strength of 3340 oersteds. The measured value 
for the sample of pyrrhotite oriented parallel to 
the field is lowest, 152 oersteds; the random 
sample is next, 220 oersteds; and the cross- 
oriented sample highest, 237 oersteds. 

These results, which are strikingly different 
from those obtained with magnetites in which 
the coercive force was independent of particle 
orientation, suggested that the coercive forces of 
random and oriented samples of magnetite as 
determined at about one-fourth saturation 
might exhibit a similar difference. Results for 
different degrees of saturation of magnetite ob- 
tained by determinations using many small 
steps in the region where the hysteresis curves 
cut the axes are shown in Fig. 5. 

The coercive force increases with the degree of 
magnetic saturation, as would be expected—14 
oersteds at Hymax.=111.5, 17.5 oersteds at 
Hiax.=198, and 19 oersteds at Hmax. =3340. 
Coercive force, however, is independent of 
particle orientation in ali cases. 

The main difference between the action for 
magnetite and that for pyrrhotite is apparently 
associated with a greater lag in demagnetization 
inherent in those particles of the latter mineral 
oriented other than parallel to the field. The 
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Fic. 5. Hysteresis curves, upper branch near origin, for 
—65+250 magnetite, different maximum magnetizations. 


position of easiest magnetization is also the 
position of easiest demagnetization for both 
minerals. 
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The results for pyrrhotite are in even more 
perfect accord than those for magnetite with the 
postulations of Sizoo® on the magnetic properties 
of single crystals in which he predicted that 
crystal orientation would affect all three proper- 
ties—induction, remanence, and coercive force. 

The results given in this paper conform to the 
modified Weber domain theory discussed in a 
previous paper of this series. The difference in 
action of orientation on the coercive force of 
magnetite and of pyrrhotite suggests that the 
study of other samples of these minerals, es- 
pecially that of normal pyrrhotite from Morro- 
Velho, Brazil, might prove very instructive. 
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Studies were made of the unequal potential distribution 
commonly observed in insulating liquids under continuous 
potential difference. The potential distribution between 
parallel plates is measured in two insulating oils for over-all 
voltages between 500 and 1500 volts, and at temperatures 
between 20 and 60°C. This distribution is produced by an 
accumulation of space charge near the electrodes which 
accumulation, however, extends measurably into the 
volume of the liquid. The amount of space charge as com- 
puted from the potential distribution increases with in- 
crease of applied voltage and decreases with increase in 
temperature. On reversal of the applied voltage, the initial 
over-all conductivity is of approximately the same value 
as that preceding reversal. The initial current on reversal 
may remain constant over an appreciable interval, say 
from 1 to 4 sec. Following the initial constant value, the 
reverse current rises sharply to a maximum from 2 to 6 


times the initial constant value, and at first rapidly and 
then slowly decreases to the same long-time values as that 
of the preceding opposite potential. In this condition, a 
reversal of the foregoing pattern of potential distribution 
takes place. 

It is suggested that two types of ion are involved in the 
space charge, one of high mobility, and the other of low 
mobility, due to the accumulation about it of an aggregate 
of neutral molecules. Oscillographic studies of the variation 
of current on reversal offer means for approximate compu- 
tations of the mobility of the ions involved. The mobility 
of the ions in a given liquid may take a wide range of values. 
As between different oils and for the same oil at different 
temperatures, mobilities vary inversely as the viscosity. 
In a low viscosity insulating oil a high value of 10-4 
cm/sec./volt/cm has been observed at 22°C. 





INTRODUCTION 


N a recent paper Whitehead and Marvin! 

have studied the nonlinear potential distribu- 
tion between the two parallel plate electrodes of 
a condenser having insulating oil as dielectric. 
They attribute the phenomenon to the accumu- 
lation of space charges in the volume of the oil 
near the electrodes. They further showed that if 
the voltage applied to the condenser is suddenly 
reversed, the reverse current passes through a 
pronounced maximum as the original potential 
pattern is built up in the opposite sense. From 
the dimensions of the condenser and the time 
required for the reversal of the potential pattern, 
they estimated that the mobility of the ions in- 
volved was of the order of 10~* cm/sec./volt/cm. 
Herzfeld,?: * using the Whitehead-Marvin curves 
and the ionization theory of gases as extended to 
liquids, arrives at the figure 10~* for the mobility. 
Two different types of ion or perhaps a wide 
range of type are thereby suggested. Schumann,‘ 
since the above work, has studied analytically 
the current variations to be expected in the 
charge and discharge currents of condensers 


1 J. B. Whitehead and R. H. Marvin, Trans. A. I. E. E- 
49, 647 (1930). 

2K. F. Herzfeld, Phys. Rev. 37, 287 (1931). 

3K. F. Herzfeld, Phys. Rev. 34, 791 (1929). 

*W. O. Schumann, Arch. Elek. (27), p. 241 (1933). 


due to the motions of layers of space charge in 
the liquid dielectric, and concludes that the 
maxima of the type reported here are to be ex- 
pected under certain conditions. This paper 
reports some further studies of the subject. 


SCOPE OF THE WORK 


Two high-grade insulating oils were examined 
with special reference to (1) the distribution and 
the variation with time, following the initial 
application of voltage, of the potential gradient 
within the dielectric; (2) the magnitudes of the 
space charges involved; (3) the relation between 
the current and the potential variations following 
the abrupt reversal of applied voltage; (4) the 
mobilities of the ions involved, with special 
reference to the reverse current maximum. 

Commercial insulating oils were used. Oil (A) 
was a clear, highly refined insulating oil, vis- 
cosity 0.23 poise at 25°C, average 30-minute 
conductivity 5X10-" mbhos/cm’*. Oil (B) was 
somewhat less highly refined, viscosity 15.0 
poises at 25°C, average 30-minute conductivity 
6X10-" mhos/cm*. The oils were received in 
sealed containers and were placed in the test 
cell without further treatment. Several other 
materials, principally waxes, were tested with 
reference to potential distribution only. The 
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results on only one of these, a commercial 
paraffin, are reported. Its melting point was 
approximately 54°C, conductivity at 68°C, 
2.6X10-" mhos/cm*. 

The principal quantities measured were: cur- 
rent vs. time, potential vs. time, and distribution 
of potential gradient in the final steady state. 
The same observations were made for an abrupt 
reversal of voltage following a foregoing steady 
state. The observations were made at each of 
three voltages, 500, 1000, and 1500 volts, and 
over the temperature range 20°C to 45°C. 


EXPERIMENTAL METHOD 


The variation of the current over the first 
few seconds following the application or reversal 
of voltage was observed with an Ejinthoven® 
string galvanometer and rotating photographic 
film, operated in conjunction with a rapid charg- 
ing, short-circuiting, discharging, and reversing 
switch, as described by Whitehead and Marvin.® 
The arrangement provides for a brief instant of 
short circuit of the condenser before the applica- 
tion of reverse voltage, for the removal of the 
geometric charge on the plates. For the values of 
the current at time intervals greater than one 
minute, a sensitive D’Arsonval galvanometer 
was used. 

The oil to be tested was placed in a carefully 
cleaned open container into which were im- 
mersed a pair of parallel brass plates, each 311 
sq. cm in area, and each completely surrounded 
in its own plane by a guard ring 3.7 cm wide. The 
back of each plate and the measuring electrodes 
were completely shielded. The insulation was 
equipped with guards maintaining leakage at a 
negligible level at all exposed points. The plate 
separation was variable up to a maximum of 
1.90 cm. 

The potential in any given plane in the oil 
between the electrodes was measured with a 
test electrode or probe. The probe potential was 
determined by a sensitive! zero method which 
avoided the passage of any charging current 
from the probe to the detecting gold leaf electro- 
scope. The test probe was a thin metallic 


®>W. Einthoven, Ann. d. Physik 12, 1054 (1903). 
6 J. B. Whitehead and R. H. Marvin, Trans. A. I. E. E. 
48, 299 (1929). 


rectangular frame 6.3 cm? wound zigzag with 
0.005 in. Therlo wire forming a grating of wide 
mesh. It was supported from a carriage moving 
in a slot attached to the electrode super-struc- 
ture. The position of the probe was read with a 
micrometer screw. The test cell was surrounded 
by an oil bath, the temperature of which was 
maintained constant by means of immersed 
electric heaters and thermostatic control. Pro- 
vision was also made for low temperature studies 
using carbon dioxide snow. The entire equipment 
including test cell and temperature bath was 
placed in a heat insulating container. 


RESULTS OF EXPERIMENTS 
Potential gradient 


In Figs. 1, 2, and 3, curves (a) indicate the 
characteristic variation with time of the poten- 
tial at a given point in the oil following the first 
application of voltage. Plate separation 1.9 cm, 
temperature 23°C, voltage 500, 1000, and 1500 
volts. These curves show that at the instant of 
application of voltage, the potential at any point 
in the dielectric is that corresponding to a uni- 
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Fic. 1. Potential time curves and potential distribution for 
500 volts applied to oil (A). 
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Fic. 2. Potential time curves and potential distribution for 
1000 volts applied to oil (A). 
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Fic. 3. Potential time curves and potential distribution for 
1500 volts applied to oil (A). 
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Fic. 4. Potential distribution for 500, 1000 and 1500 volts 
applied to oil (B). 
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form gradient between electrodes. The potential 


’ then rises (or falls) rapidly from this value and 


arrives at a steady value only after a period of 
twenty minutes or more. The shape of these 
curves and the time intervals involved suggest 
an accumulation of charges of opposite sign in 
the neighborhood of each electrode, that the 
process takes place slowly, and that very slug- 
gish ions of low mobility are involved. 

Curves (b) of Figs. 1, 2, and 3 show the final 
potential distribution in the steady state for each 
of three values of applied voltage. The curves of 
Fig. 4 show the potential distribution for oil B. 
The curves indicate a symmetrical variation 
about the central plane with relatively high 
values of gradient in the neighborhood of the 
electrodes. This electrode gradient is greater the 
higher the applied voltage. Obviously the gradi- 
ent at the center of the dielectric is correspond- 
ingly less than the initial uniform value. The 
shapes of these curves also suggest volumetri- 
cally and symmetrically distributed space charges 
near the electrodes and with opposite signs. 











SPACE CHARGE IN 

Approximate values of potential gradients 
may be taken directly from the curves, and 
may be used for a computation of the total 
amount of space charge accumulating on either 
side of the central plane. 

If F, and F, are the measured gradients, in 
electrostatic units, at one electrode, and in the 
center of the dielectric respectively, then by 
Poisson’s equation, we have for the total space 
charge per cm? of plate area between plate and 
central plane of the dielectric: 


af K 
a= f pdx —on F,,) statcoulombs/sq. cm 
a T 


or, since F, and F, are measured in volts, then 
the charge expressed in coulombs is: 


K 1 
Q=—(F— P)( ) 
4n 9x 10" 


coulombs 





sq. cm of plate 


Values as measured and computed in this way 
and corresponding to the data of Figs. 1, 2, and 3 
are given in Table I. The values are subject to 
an error of +25 percent because of the uncertain- 
ty in the evaluation of the potential gradient. The 
final values of the space charge are based on the 
mean of the values of the gradient at the two 
electrodes. It is seen that the total space charge 
is relatively small. However, it increases with 
the applied voltage, and, as the over-all average 
gradient is relatively low in these experiments, at 
higher values of gradient much greater values of 
the space charge are to be expected. 

If in the steady state of the dielectric, the 
applied voltage is removed, the potential of the 
test probe falls more or less slowly. The rate of 
fall is lower the higher the applied voltage and 


TABLE I. Space charge. Oils A and B. 








GRADIENT (volts/cm) SPACE CHARGE/cm? 





SPECI- Stat- 
MEN Vo.tts Temp. , Ave. Faq F, Coulombs coulombs 
Oil (A), Plate Separation 1.9 cm 
A 500 22°C 1000 130 1.55x10-" 0.465 
B 1000 = 22.5 1420 180 2.2 «107 .662 
C 1500 23 3300 380 5.1 «107 1.55 
D 1500 48 2675 510 4.0 «10° 1.22 
Oil (B), Plate Separation 1.1 cm 
E 500 20 2450 180 3.8 x10 1.15 
F 1000 20 4450 300 7.3 x10-" 2.21 
G 1500 20 4475 550 7 x10" 2.1 
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the lower the temperature. The fall is more rapid 
if the geometric charge is removed by a brief 
instant of short circuit. This behavior is in 
general accord with the behavior indicated by 
the numerical data, and on the assumption that 
the space charges are made up of mobile ions. 

If in the steady state under applied voltage, 
the latter is removed, the electrodes are short- 
circuited for a brief instant, and then connected 
through a sensitive galvanometer, no current is 
observed. In other words, there is no discharge 
current corresponding to the slowly decaying 
charging current commonly observed. This type 
of behavior has been frequently noted. It indi- 
cates that the accumulation of space charges 
near the electrodes is in large measure restrained 
from recombination by corresponding bound 
charges on the electrodes. The slow return’ to the 
neutral condition through diffusion and slow 
recombination is by a low value conduction 
current in the external circuit as the opposite 
space charges slowly diffuse together. These two 
are in opposite sense and their difference is so 
small as to be measured only with high sensitivity 
equipment. 

The sign and volumetric distribution of the 
space charges was shown by freezing the oil 
while under stress, with carbon dioxide snow. 
The frozen oil was then removed from between 
the electrodes cut in half in a plane normal to the 
field, and each half placed, successively in a 
Faraday pail, giving the sign of the charge and 
approximate value of potential. The experiment 
is difficult as the oil quickly loses both solid 
state and stored charge. However, the results 
indicated clearly substantial charges in each of 
the halves, corresponding to the potential dis- 
tribution curves of Figs. 1, 2, 3, and 4. A volumet- 
ric distribution of charge rather than a surface 
layer at the electrodes is thereby indicated. 

The potential distribution in paraffin was also 
measured (see Fig. 5). The paraffin was melted 
into the condenser, the voltage applied while 
it was in the liquid state until the steady condi- 
tion was reached. It was then allowed to solidify, 
voltage being maintained on the electrodes. 
After solidification, the potential distribution 


7J. B. Whitehead and S. H. Shevki, ““D-C Clean-up 
in Insulating Oils,” presented Summer Convention, ~ 
A. I. E. E., June 24-28, 1935. 
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Fic. 5. Potential distribution for paraffin. 1500 volts 
applied. 


as shown in Fig. 5 was measured. The unsym- 
metrical distribution indicates that the space 
charge is practically entirely of one sign and 
suggests that the ions of one sign pass freely to 
the electrode whereas those of the other sign do 
not. Following these observations, the solid 
block of dielectric was removed from the con- 
denser and tests with the Faraday pail indicated 
a negative charge corresponding to the location 
and shape of the potential distribution curve. 
The paraffin was also cut into two pieces in a 
plane parallel to the plates and the sum of the 
potentials observed for each half was equal to 
that of the whole piece before cutting. Experi- 
ments of this character must be performed 





WHITEHEAD AND E. E. 


MINOR 





quickly as the paraffin loses its stored charge 
rapidly unless carefully insulated. These observa- 
tions also indicate clearly a distribution of the 
total space charge through the whole volume of 
the dielectric. 


Behavior on reverse potential 


If in the steady state under applied voltage, the 
voltage is suddenly removed, the plates short 
circuited for a brief instant, and then the same 
value of potential is applied in the opposite 
sense, the variation of the current with time 
may be followed by photographic means, using 
the string galvanometer or the oscillograph.* 
Photographs of the reverse current show that 
the current always remains constant for a few 
seconds after reversal, then rises quickly to a 
peak value and then falls away at much the 
same rate, or sometimes very slowly, to a final 
steady value. A tracing from a typical photo- 
graphic exposure is shown in Fig. 6. The reverse 
current curves for oil (A) shown in Figs. 7 and 8 
are plotted from the photographs; the upper 
curves show the corresponding potential varia- 
tion at a selected point in the dielectric. The 
curves indicate that the potential on reversal 
passes through practically its entire variation or 
has arrived at its full value when the sharp 
maximum in the current has passed, and that it 
varies most rapidly while the current is passing 
through its most rapid variation. Numerical 
data taken from the curves of Figs. 7 and 8 are 
given in Table II. They indicate that following 
reversal, the current remains constant in value 
for a second or more. This value is usually from 
1 to 2 times the value of the long time steady- 
state current before reversal of voltage, but never 


8 J. B. Whitehead, Trans. Electrochem. Soc. 65, 83 (1934). 
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TABLE II. Jonic mobilities. Oils (A) and (B). 
CURRENT (amp. X 10°!°) Time (Sec. Tora SPACE 
TEMP.. before after reversal current max.) CHARGE MOBILITY 
SPECIMEN VoLts (°C) reversal initial max. start crest (coulombs) (cm/sec./volts/cm) 
Oil (A), Plate Separation 1.9 cm 
B 1000 22.5 4.9 5.0 20 14.0 31.0 6.8 x 1078 1.06 x 10-4 
c 1500 23 7.0 9.0 30 4.5 15.0 15.8107 1.04 10™4 
B’ 1000 44 10.0 11.5 23 2.6 11.0 5.21078 2.15 x10~4 
ity 1500 47 13.0 12.5 45 1.1 6.1 11.7107 2.4 x10~4 
Oil (B), Plate Separation 1.1 cm 
H 1500 58 102.0 106.0 189.0 35 600 20.7 x 10-8 2.1 «K10~° 
H’ 1500 58 96.5 63.0 130.0 40 600 19.5 x 1075 2.1 «10°° 








greater. The reverse current maximum is from 
2 to 6 times that of the earlier steady-state cur- 
rent and following the sharp maximum may 
decrease slowly over a long time to the final 
steady value: The maximum increases with the 
potential gradient and with temperature. The 
reverse Maximum is greatest in magnitude when 
observed after the space charge has reached 
equilibrium condition. The higher the tempera- 
ture and the higher the applied voltage, the 
shorter the duration of the initial constant cur- 
rent and the sharper the current maximum on 
reversal. For applied voltages less than 500 volts, 
the reverse maximum is very low and detected 
only with difficulty. Values of space charge com- 
puted by integrating the area of the reverse 
current maximum, i.e., the area of above the 
long time steady (conduction) current, are 
of the same order of magnitude as those com- 
puted from the potential distribution. 


Computation of ionic mobility 


The potential distribution curves reverse 
themselves accurately on reversal of applied 
voltage. The time required for such reversal is 
clearly indicated by the reverse current and 
potential time curves of Figs. 7 and 8. Knowing 
the distance between electrodes, and assuming 
that the new potential distribution on reversal 
is caused by the passage of the space charges to 
the opposite electrodes, we have a simple over- 
all means of determining the average mobility 
of the ions making up the space charges. Ac- 
cordingly we have: k=d/E-t, in which k=mo- 
bility in cm/sec./volt/cm, E=gradient in volts/ 
cm, d=average length of path of ions in cm, 
t=time in sec. for reversal. 


In applying this formula, the most uncertain 
factor is ¢ for in some cases the approach to the 
final steady value after reversal is very slow; 
however, the steep slopes of the current maxima 
indicate that the greater part of the space charge 
must change position more or less rapidly. A 
more accurate determination appears probable 
in using for the elapsed time, the interval be- 
tween the initial rise of the current maximum 
and the maximum itself, and assuming that at 
the latter instant the space charges have trav- 
ersed half the distance between the plates. 
Values of mobility so computed are given in 
Table II. It is seen that within the range meas- 
ured, they are independent of the potential 
gradient, increase with temperature, and are 
roughly 100 times greater in oil (A) than in 
oil (B), corresponding to the large difference in 
viscosity of the two oils. It is to be noted also 
that the order of magnitude of the mobility in 
oil (A) is 10-* cm/sec./volt/cm, corresponding 
to the earlier measurements of Whitehead and 
Marvin. It may also be noted that at the maxi- 
mum of reverse current, the reverse potential 
distribution has been largely established. This 
means that at this time a portion of the oil at 
least is subjected to the counter-gradient of the 
polarization caused by the space charge. This 
means that the resultant gradient is less than 
that corresponding to the applied potential and 
the average mobilities computed above may have 
somewhat higher values. 

The mobility may also be computed from the 
value of the reverse current maximum and the 
total volume of space charge, again assuming 
that in reversal the total space charge reverses 
its position. Two values of the space charge are 
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Fic. 8. Variation of potential and current on reversal at 1500 volts. 








| 
| 
4. 

















eet 


| 
| 
| 
| 





MINOR 


| 
+ 
AT 0461cM 








1S00V, — 22.5°C 
1500V, — 47°C 














| 


(ee Neen tee 


j—________} 








oe 


oe mae 





— 
° 
x 


E 
PROBE AT 0.461 CM 


x lwoov, ~— 23°C 
O 1wooVv,—44°C 


oe 
. en 





— 
| 
| 


— PoTentiac 
| 


——-——1 


a 


AND E. 
--~-~-~~1 
Time, SECONDS 


| 
80 








PotenTiar 





























— 
—_ aa Ne 
| 
| 
~ 





—— 
| 


~ 
—_ 
< 
ea) 
_ 
—_— 
i) 
‘ 

—_ 
— 
_ 
_ 
> 
> 





B. 











| 
iy | 
a 
| | 


Fic. 7. Variation of potential and current on reversal at 1000 volts. 
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available, one as taken from the potential dis- 
tribution curves and the other by integrating 
the area which lies above the steady state 
(conduction) current. The two values of space 
charge so computed are of the same order of 
magnitude, but do not agree closely. The mobil- 
ities so computed for oil (A) are about one-tenth 
the values computed by the first method and as 
given in Table II. A possible explanation is the 
recombination of the ions of opposite sign, as 
the two space charges must pass each other in 
moving to the opposite plates. 


DISCUSSION 


Two features of the foregoing observations 
stand out as not being in accord with present 
theories of conduction in liquid dielectrics. They 
are: (1) the reverse current maximum, and (2) 
the initial relatively low constant value of 
current on reversal. 

When in the steady state, the current voltage 
relations in a pure insulating liquid are very 
similar to those in gases, namely, initial linear 
relation, followed by saturation, and then rapidly 
rising current due to secondary ionization. 
These and other relations indicate strongly that it 
should be possible to extend the ionization theory 
as developed for gases also to the phenomena of 
liquids.* Unequal potential distribution and 
space charges are found in gases and are in ac- 
cord with gas theory. The space charges are due 
to space variations in ionic density, but the ions 
themselves are in uniform motion towards the 
electrodes. When the applied voltage on parallel 
plates is reversed, the current after reversal, in 
accordance with gas theory, should be approxi- 
mately of the same value as before reversal, 
and should not vary with time. Consequently, 
the pronounced maxima found in the reverse 
current behavior of insulating liquids is ap- 
parently not accounted for by the laws explain- 
ing space charge and current in gases. The initial 
constant current on reversal is apparently in 
accord with the behavior to be expected from 
gas-laws. Its order of magnitude is the same as 
the preceding long-time steady state current. 
Usually it is slightly greater, and this is readily 


*G. Mie, Ann. d. Physik 26, 597 (1908). 
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explained by the fact that immediately after 
reversal, the space charge polarization and the 
applied voltage are now both in the same direc- 
tion so that the over-all gradient at the beginning 
of reversal is somewhat higher than that cor- 
responding to the applied voltage. On the other 
hand, this initial current, although constant over 
an appreciable interval, then merges into the 
rapidly rising reverse current maximum and 
there appears no obvious reason why the two 
should be separated in seeking for an explanation. 

In the neutral state, even when highly purified, 
insulating liquids always show some conduc- 
tivity albeit of very low value, due to traces of 
impurities, radioactive influence, etc. The total 
number of ions present rises to equilibrium condi- 
tions as related to recombination and diffusion. 
Some of the free ions attract to themselves 
clouds of neutral molecules, thus constituting 
large and sluggish ions. Under normal electric 
gradients, an ion may carry with it as it moves a 
part of this surrounding cloud. As the gradient 
increases, however, it may be expected finally to 
shed the cloud completely. 

When a potential difference is applied, if not 
too high, there will be motion of both the old 
heavy ions and newly generated free ions towards 
the two electrodes. The old heavy ions will move 
much more slowly and if they are sufficient in 
number a practically steady initial value of cur- 
rent may result as often observed.*: '® As the 
ions arrive at the electrodes, some of them are 
discharged, but others of them are not dis- 
charged, retaining enough of their neutral 
envelopes to prevent complete arrival at the 
electrode surface. Space charges will, therefore, 
accumulate near the electrodes due partly to 
the steady withdrawal of freshly formed ions 
from the central regions in accordance with gas 
behavior, and also due to the accumulated and 
more or less stationary heavy ions with their 
surrounding clouds of neutral molecules. As 
the space charge increases and the internal 
gradient decreases, there is a greater and 
greater tendency for an increase in the number 
of these heavy ions with relatively wide distribu- 
tion through the volume of the liquid. Both of 
these elements of space charge are necessary for a 


1” J. B. Whitehead, Trans. A. I. E. E. 50, 692 (1931). 
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complete explanation of the phenomena as 
observed. 

The long-time steady state current is due 
principally to the constant supply of newly 
formed more mobile ions. The initial reverse 
current is due only to the newly generated 
mobile ions, the conductivity for this brief 
period being the same as the final conductivity of 
the preceding opposite period. The reverse 
current maximum is apparently due _prin- 
cipally to the space charge stored up in the 
large sluggish ions, consisting of molecular ag- 
gregates. No such maximum is to be expected 
from the new more mobile ions. A difficulty in 
this view exists, however, in the time interval 
which elapses before the rise to the maximum 
begins. Considered merely as mobile ions, we 
would expect them to take up their new direc- 
tion and velocity immediately on reversal of 
the field. A possible explanation suggests itself 
as follows: The heavy ions with their envelopes 
of neutral ions pack more or less tightly in the 
region nearest the electrodes. Because of their 
size, the mobilities of these aggregates must be 
very low at the start, adding practically nothing 
to the conduction current. However, under the 
continued action of the field, particularly in the 
central region where the density of the heavy 
ions is less, and under the added influence of the 
counter-polarization, the ions, as they get into 
motion, begin to shed their neutral envelopes, 
thereby assuming a steadily increasing velocity. 
This process communicates itself in turn to the 
aggregates lying nearer the plates until finally 
the whole volume of heavy ions have shed their 
neutral molecules and assume their maximum 
possible mobility, thus setting up the reverse 
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current maximum as observed. The accumula- 
tion of the opposite space charge then causes a 
decrease in velocity and a further accumulation 
of neutral envelopes and a return to the heavy 
sluggish state. 

In connection with the value 10~* cm/sec./ 
volt/cm for the mobility in oil (A), it is to be 
noted that this figure is of the order of that cor- 
responding to the mobilities of ions in electro- 
lytic solutions. This is in accord with theoretical 
discussions by Gemant," Nikuradse,” and others 
who ascribe an electrolytic character to the 
conductivity of insulating liquids. Nikuradse 
in fact has mentioned even higher mobilities of 
the order 10-* for certain insulating liquids. 
Hofmann," using an optical method and cinema 
camera, has shown recently that the mobility of 
the ions in a transformer oil immediately fol- 
lowing the first application of voltage is of the 
order 10-* cm/sec./volt/cm. All these figures 
strongly suggest that the ions involved in the 
conductivity of insulating liquids in certain 
cases at least are electrolytic in their origin. 
However, it is to be noted that for the oils 
reported in this paper, as well as the others 
referred to, all have substantially higher values 
of viscosity than those for the usual electrolytic 
solutions, which would call for a lower mobility 
if the ion is of simple electrolytic character. It 
appears probable therefore that in the case of 
insulating liquids, other elements may be present 
than those included in the common theory of 
electrolytic dissociation and conduction. 

"A. Gemant, Elektrophysik die Isolierstoffe (Julius 
Springer, Berlin, 1930). 

12 A. Nikuradse, Das fliissige Dielektrikum (J. Springer, 


Berlin). 
183 R. Hofmann, Zeits. f. Physik 92, 739 (1934). 
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Dynamical similarity laws of the high pressure mercury 
vapor discharge have been set up. These laws are based on 
principles which determine characteristics of a high pres- 
sure discharge in thermal equilibrium. Characteristics of 
discharges can be determined by measurements on models 
if the model contains the same mass of mercury gas per cm 


and is supplied with the same wattage per cm of length of 
tube as the main construction. The temperatures of the 
walls of discharge tube should be the same. In formulating 
the above requirements self-absorption has been neglected. 
Experimental tests have shown the validity of the similarity 
theorems. 





§1. INTRODUCTION 


REQUENTLY it occurs that the engineer, 

though knowing the physical principles which 
govern the processes in a newly constructed 
device, is unable to predetermine its character- 
istics by calculation because of the mathematical 
difficulties involved in this problem. 

In such circumstances he has recourse to 
experiment, determining the interesting features 
of his new design by means of models. In order 
to apply the data yielded by model research to 
the main construction C the model M should 
correspond with C in all actual details. It is 
often impossible, however, to make the dimen- 
sions of C and M equal. Thus there arises the 
problem, how to effect a compromise which will 
ensure, that the processes of C and M are 
analogous in such a way that measurements 
obtained from M are conclusive with respect to C. 

A general method of treating this problem may 
be outlined as follows: If C and M are similar, 
it should be possible to choose such different 
units, measuring corresponding quantities in C 
and M, that the differential equations and 
boundary conditions determining the process 
investigated agree numerically in both devices. 
Then, the results of investigation will also agree 
numerically and the conversion factor follows 
from the dimension of the measured quantity. 
The characteristic quantities of zero dimension, 
in particular, should have the same value in C 
and M. 


§2. THE PHYSICAL PRINCIPLES OF THE MERCURY 
HIGH PRESSURE DISCHARGE 


From the preceding section it is apparent 
that similarity considerations are possible only 


if the physical principles of the process are 
known. It is necessary therefore, first to give a 
short treatment of the physical theory. The 
principles determining the characteristics of the 
high pressure discharge and the similarity laws 
obtained therefrom are known from compre- 
hensive investigations of Elenbaas.'~* 

Contrary to the low pressure discharge the 
high pressure discharge is in thermal equilibrium. 
The temperature of the gas determines the 
velocity distribution of atoms, ions and electrons; 
it determines the degree of ionization and by this 
the electric conductivity of the positive column 
and thus the wattage supplied by a given cur- 
rent. It further determines the wattage dissipated 
by radiation, heat conduction and heat convec- 
tion. From the differential equation of the tem- 
perature, deduced in the following two sections, 
the similarity laws of the positive column will 
be readily found. We confine ourselves to con- 
sidering 1 cm of an infinitely long column. The 
following notation will be used: 


Notation 


r,% 


Significance 

Cylinder coordinates 
Inner radius of the tube 
r/o 
Time 
Absolute temperature 
Pressure 
Electric field in the direction of 

the axis (gradient) 
Current density 
Density of the mercury vapor 
Hg mass in 1 cm tube length 
Thermal conductivity 
Mobility of electrons 


wa Foe Q Where 


1 W. Elenbaas, Physica 1, 211 (1934). 
2 W. Elenbaas, Physica 1, 673 (1934). 
3 W. Elenbaas, Physica 2, 45 (1935). 

4W. Elenbaas, Physica 2, 169 (1935). 
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Degree of ionization a 
Ionization energy of the mercury 
atom €;=1.65-10 “erg 
=k-120,000°K 


Charge of 1 g of Hg ions q= 1.445-10" e.s.u. 
Boltzmann constant k=1.37-10-"* erg/grad 
Constant of Saha’sformula K=0.335 g cm™ sec™? grad~*? 
Constant of Langevin’s formula Cc 


The notation given above refers to M; quan- 
tities referring to C are included in brackets. 


(a) Ionization and current density 

Saha’s® formula determines the degree of 
ionization: 

oe P/(1—a?) = K-T*/e-«i/k7, (1) 

For all temperatures attainable in practice 
a is <1, because ¢«;/kT will remain large even 
at the highest temperatures. 

Using the laws of an ideal gas we assume P to 
be proportional to m/r,?. Neglecting a? in the 
denominator of (1) there follows: 


oem/re= S(T), (2) 


whereby S(T) is a function of the temperature 
only. 

From the mobility B of the electrons the 
current density: 


i=a-p-q-B-G (3) 


is obtained, neglecting the mobility of ions. 
B is inversely proportional to the concentra- 
tion p. According to Langevin® we may assume: 


B-p=C/T'. (4) 
Combining (4) and (2) we obtain from (3): 
i= (Gro/m')-C-qgLS(T)/T]}. (5) 


(b) Energy balance and differential equation of 
the temperature 


Consider a cylindrical sheet of radius r=r9-x 
and thickness dr=rodx. The wattage supplied is 
according to (5) 


dL=i-G2rrdr 
= (G’ro3/m')[C-q(S/T)!2axxdx]. (6) 


The same wattage must also be dissipated. 


5M. N. Saha, Phil. Mag. 40, 472 (1920). 
® Langevin, Ann. d. chim. et phys. 81, 


5, 245 (1905). 
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Neglecting the contribution of convection, only 
thermal conduction and radiation will be taken 
into account. The wattage given off by thermal 
conduction amounts to: 


AWi=(aera(T)—)ar="(2220—-)ax (7) 


The radiation, when dealing with the phe- 
nomenon of self-absorption, is subject to very 
complicated laws. Fortunately the greater part 
of energy radiation is coming from transitions 
between excited states of the atom. As mercury 
atoms in those excited states exist only in very 
small concentrations, we may assume that self- 
absorption is of little influence. Then the radiant 
energy will be proportional to the number of 
excited mercury atoms: 


dW2=m2rxdxF(T). (8) 


In this formula F(T) is a function of the 
temperature only which actually contains the 
excitation energy ¢, in the Boltzmann factor 
e~«/kT, Combining (6), (7), (8) we obtain from 
the energy balance: 


dL=dW,+dwW, 


the differential equation: 


d dT 
© (xo(T)— +mF(T)x 


dx 
_ Gr! S( T) 
(co) +) 
m} 

This differential equation determines the tem- 
perature as a function of x and of the parameters 
m, G, fo. Knowing the temperature distribution, 
we may calculate from (5) the current density, 
from (6) the wattage supplied, from (7) the heat, 
transported by conduction, from (8) the energy 


radiated. We see that the characteristic prop- 
erties of the column are entirely determined. 





$3. THE SIMILARITY PRINCIPLES OF THE 
Mercury HIGH PRESSURE DISCHARGE 


As we cannot solve Eq. (9), the characteristics 
of a model should be experimentally determined. 
It is the purpose of similarity dynamics to 
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Fic. 1. Product gradient-(diameter)? as a function of 
the power input, measured with tubes of different diameters 
d. The similarity laws require this product to be inde- 
pendent of d. 


predetermine from model research the particulars 
of the main construction with other parameters 
(m), (G) and (ro). It has been indicated in the 
introduction how this is effected. The parameters 
of M should be chosen in such a manner, that 
Eq. (9) agrees numerically for M and C. It is 
apparent, that this only holds, when: 


m=(m) and G?r.3/m!=(G)?(ro)3/(m)'. (10) 
Then we obtain from the Eqs. (6), (7), (8): 
L=(L), Wi=(Wi), We=(We). (11) 


Hence we may state the theorem: 

High pressure discharges are similar if they 
agree in the quantity m of mercury gas and in 
the wattage supplied both per cm of length of 
the tube. In similar discharges the squares of the 
gradients are inversely proportional to the 
third power of the tube radii. Similar discharges 
dissipate the same wattage (per cm of length 
of the tube) by conduction as well as by radiation. 

To the two requirements for similarity, 
mentioned above, should be added that the 
boundary conditions of (9) i.e., the temperature 
of the bulb, shall agree. To guarantee this, for 
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two discharges of different diameter, the smaller 
one should be cooled more thoroughly. For- 
tunately, the gradient of the column has proved 
to. be practically independent of the external 
temperature. (When heating up to 300°C 
Elenbaas‘ found an increase of approximately 
three percent.) The influence of external tem- 
perature thus will be neglected in the following 
considerations. 

To test our theorem of similarity it should be 
proved that G*r,? depends on L and m but not on 
ro. This similarity law is well confirmed by ex- 
periment as diagrammatically shown in Fig. 1. 
In the range of diameters between 8 mm and 39 
mm the quantity G-7,*/? appears to be the same 
function of Z and the parameter m within the 
accuracy of measurements. 


§4. GENERALIZATION OF THE CONCEPTION OF 
SIMILARITY 


Integrating Eqs. (7), (8), (9) over x we 
obtain: 


Wi=A, a 
W.=B-m, b?+(12) 
L=A+Bm = (G?ro3/m!)-D. Cc 


Here A, B, D are integrals, depending on the 
temperature distribution only. They will be 
considered as functions of the “effective average 
temperature Ts, of the path of discharge, as- 
suming by hypothesis, that the value of these 
integrals at a given effective temperature does 
not actually depend on other parameters, in 
particular on m. Considering the agreement of 
effective temperatures as similarity in a broader 
sense, the following similarity laws are obtained. 





L/(A+Bm) =(L)/[A+B(m)], (13) 
G*rom-* (G)?(r9)3(m)-3 
= ; (14) 
A+Bm A+B(m) 
Thus L/G?ro'm-* = (L)/(G)*(r0)3(m)~! (15) 
and further: W,=(W)), (16) 
W2/m=(W2)/(m). (17) 


Utilizing these equations we may state 
the generalized theorem: High pressure dis- 
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Fic. 2. G*d*/m' as a function of m. Crosses: measure- 
ments with similar discharges, taken from the curves of 
Fig. 1. Drawn line: the linear function A’+B’m with: 
A’=8.5-10-4, B’=5.75-10-*. 


charges are similar if they agree in the quantity 
L/G*rem-'. 

In similar discharges the gradient is deter- 

mined by 

G?(2ro)*/m'=8(A+Bm)/D=A’'+B’m. (18) 
Similar discharges dissipate the same wattage by 
conduction whilst the radiation will be propor- 
tional to the (different) values of m. 

This theorem may be tested simultaneously 
with the determination of A’ and B’ by repre- 
senting G?(2r,)*/m} as a function of m for similar 
discharges. Fig. 2 drawn from the same measure- 
ments as Fig. 1 indicates that Eq. (18) is fully 
confirmed. The drawn line corresponds with: 

A’=8.5-10', B’=5.75-10', D=6.25. 
These values have been proved to hold between 
the limits: 


20 watt/cem < L <50 watt/cm, 


? mg/cm <m<12 mg/cm. 


$5. APPLICATION OF THE SIMILARITY THEOREM 
TO LiGHt EMIssION 


Theory 


Eqs. (13), (16), (17) yield the following re- 
markable conclusions: 

(a) From the wattage supplied and radiated, 
an efficiency factor with respect to the radiation 
may be defined by 


n=energy radiated /energy supplied. 


According to (12) a, b, c one obtains 


n= B'm/(A'+B'm); (19) 
n is a monotonous increasing function of m. As 
far as the similarity laws hold exactly, the 
spectral energy distribution of the radiation 
should agree for similar discharges. Thus the 
“economy”’ defined as the intensity (in a certain 
spectral frequency) divided by the wattage 
input should be proportional to 7. 

(b) Consider similar discharges of equal power 
input (L)=L. Then (m)=m, as is apparent 
from Eq. (13). For such discharges (W)=W 
by Eq. (17). For similar discharges of equal 
power input per cm, the energy radiated per 
cm is the same. The surface brightness of such 
discharges will be inversely proportional to the 
radius of the tube. Insofar as the laws of similar- 
ity hold, there is no limit to the extent to which 
the surface brightness may be increased by 
decreasing the radius of the tube. 


Experiment 


(a) Fig. 3 illustrates the measured economy 
compared with those calculated according to 
(19). Qualitatively the curves agree well. The 
measured lack of efficiency, increasing with rising 
values of m may be understood from the neglected 
influence of self-absorption. Because of the 
influence of self-absorption the spectral energy 
distribution of similar discharges may differ 
considerably. Interesting features concerning the 
spectrum in the case of extremely high pressures 
have been investigated by De Groot.’ 

(b) The luminous output of similar discharges 
should be independent of 7». Fig. 4 shows that 
this conclusion is experimentally well confirmed. 


7W. de Groot, De Ingenieur 50, E92 (1935). 
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Fic. 3. Economy E (arbitrary units) of the green 
mercury line \=5461A as a function of the power input. 
The units are chosen in such a way that for m=0.75 mg/cm 
[L=(8/D)(A’+B’m)=10 watts/em] the economy E 
=const. 7=1 both theoretically and experimentally. The 
difference between the measured curve and the calculated 
function may be ascribed to self-absorption. 


The slight decrease of intensity at small radii 
is again the result of self-absorption. 

If a high efficiency, combined with high sur- 
face brightness is desired, the preceding con- 
siderations direct us (a) to large values of m, 
(b) to small radii 79; both arguments point to 
extremely high pressures. Proceeding in this 
direction C. Bol designed discharge tubes of 3 mm 
diameter, bearing a pressure up to 300 atmos. 
Operating these tubes with a power of 600 
watt/cm van Alphen measured a luminous flux 
of 30,000 lumen/cm. 

To prove the applicability of our similarity 
theorems to such extreme constructions, we will 
compute the radiation according to (19), using 
as model a discharge of: 


ro=1 cm, L=40 watt/cm, m=3.1 mg Hg/cm. 


According to Elenbaas!:? this model has a 
vapor pressure of accurately 1 atmos. and a 
light output: 
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Fic. 4. Intensity of the green mercury line (arbitrary 
units) for similar discharges (m= 12 mg/cm of tube length) 


as a function of the diameter. The slight increase may be 
ascribed to self-absorption. 


W = 1650 lumen/cm. 


Hence the light output of the main construction 
will be: 





(L) (n) (L) B’(m) A’+B'm 
OU, cers uncicgal Wi Senenavs : : 
B'm 


n L A'+B’'(m) 
Introducing (L) = 600watt/cm,m = 30mg,Hg/cm, 
one obtains: 


(W) = 34,200 lumen/cm, 


which is approximately 14 percent higher than 
the measured output. Estimating the diameter 
of the contracted discharge to approximately 
1 mm, v. Alphen computed an average brightness 
of 30,000 candles/cm?. For comparison the sur- 
face brightness of some other sources of light 
is listed below. 


Tungsten filament of an incandes- 
cent lamp 200— 3000 candles/cm? 

Crater of carbon arc 9000-17,000  ‘* /cm? 

High intensity arc 40,000 “ /em? 

Sun, as observed from the surface 
of the earth 


165,000 * /em? 


§6. VALIDITY OF THE SIMILARITY THEOREM 


In the experiments of Elenbaas the tube 
diameter varied between d=8 mm andd=39 mm, 
the mass of mercury varied between m= 3 mg/cm 
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and m=12 mg/cm. The extreme cases, investi- 
gated by Elenbaas‘ and Bol,® are: 


Elenbaas: 

d=39 mm, m= } mg/cm p=0.045 atmos. (under operating 
cond.) 

d=8 mm, 
cond.) 


m=3 mg/cm p=4.2 atmos. (under operating 


Bol: 
d=3 mm, m= 
cond.) 


30 mg/cm ~=300 atmos. (under operating 


The operating pressures of Elenbaas are spread 
in a range of about one hundred-fold. At pres- 
sures smaller than a few cm Hg no thermal 
equilibrium between atoms and electrons will 
exist. The discharge gradually merges into the 
low pressure type, for which other similarity 
principles have been found. On the other hand 
however, the upper limit does not appear to be 
reached with 4 atmos. 

Two further processes, neglected until now, 
will become more and more effective at higher 


concentrations of mercury. 
I 


(a) The self-absorptisa——_____ 

Figs. 3 and 4 show, that self-absorption plays 
a noticeable role in all high pressure discharges. 
Hence the intensity of radiation can only ap- 
proximately be computed by means of similarity 
laws. An actual improvement of the theory in 
this direction seems to be impossible as the 
laws governing self-absorption are extraordi- 
narily complicated. Experiment shows that the 
electrical characteristics of the column, es- 
pecially Eq. (18) for the gradient are hardly 
influenced by self-absorption. 


(2) The heat convection 


From the energy dissipating processes, only 
heat conduction (independent of m) and radia- 


5 C. Bol, De Ingenieur 50, E91 (1935). 
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tion (proportional to m) have been taken into 
account. Both processes yield a wattage of long 
tubes proportional to their length. Estimating 
the influence of convection, an additional term 
is found, proportional to m? and independent of 
the length of the tube. This is sufficient reason 
to neglect its influence on an infinitely long tube. 
It is however not at all improbable that this 
effect becomes noticeable on the tubes made in 
practice. In that case Eq. (12c) should read: 


A+Bm-+ Cm? = (G*r,3/m')D 


and for the gradients the following formula 
would hold 


G=m'(A'+B’m+C'm?)}/r9}. 


The measurements of Elenbaas do not exclude 
the existence of such a quadratic term but 
certainly it is not of great importance. 

Apart from quantitative agreement with 
respect to the radiation the similarity laws for 
the gradient have proved to be correct even for 
the super high pressure discharge of Bol. Prob- 
ably they will hold in the entire domain of pres- 
sure from a few cm up to 300 atmos. 

(A brief bibliography is appended.) 


BIBLIOGRAPHY 


Similarity laws of the mercury low pressure discharge. 
R. Holm, Wiss. Veréff. Siemens K. 3, 159 (1923). 
Thermal ionization of arcs. 
L. S. Ornstein, Physik. Zeits. 32, 517 (1931). 
L. S. Ornstein, H. Brinkmann, and A. Beunes, Zeits. f. 
Physik 77, 72 (1933). 
A. v. Engel and M. Steenbeck, Wiss. Veréff. Siemens K. 
10, 2, 155 (1931). 
A. v. Engel and M. Steenbeck, Wiss. Veréff. Siemens K. 
12, 1, 76 (1933). 
A. v. Engel and M. Steenbeck, Wiss. Veréff. Siemens K. 
A. v. Engel and M. Steenbeck, ‘Naturwiss. 19, 212 (1931). 
R. Mannkopf, Zeits. f. Physik 76, 39 (1932). 
R. Mannkopf, Zeits. f. Physik 86, 161, (1933). 
H. Witte, Zeits. f. Physik 88, 415 (1934). 
Distribution of energy in the radiation of the mercury column. 
H. 1933). and M. Pirani, Zeits. f. tech. Physik 14, 393 
(1933). 





DECEMBER, 1935 


PHYSICS 


VOLUME 6 
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A method is described for the construction of electrical 
models with which to solve problems of the flow of liquids 
through porous media under the action of gravity. With 
these models the shape of the free surface and the extent 
of the surface of seepage are determined simultaneously 
with the potential distribution within the flow system. 
Four examples are given of the application of the method 
to the problem of the seepage of water through dams. In 


two cases the faces of the dams were assumed to be sloping, 
and in the others they were taken as vertical, the free sur- 
face and velocity distributions at the faces of one of the 
latter cases checking closely with those derived from an 
analytical solution of the problem. The method of applying 
these models to the study of more complex gravity flow 
systems is briefly indicated. 





INTRODUCTION 


ROBLEMS involving the flow of homo- 

geneous fluids through porous media under 
the action of gravity in general require very 
difficult analytical methods for solution.' Fur- 
thermore the problems may become entirely 
intractable mathematically when the geometry 
of the flow system takes on only a reasonable 
complexity such as is involved in many systems 
of practical engineering interest. For this reason 
small scale models duplicating the original in 
their physical make-up have been used in the 
study of such problems. In view of the difficulties 
involved in making such scale models as well as 
the inherent limitations in the detail and accu- 
racy available, we are presenting here a method 
based on the electrical analogy of the flow of 
fluids through porous media which is simple but 
nevertheless capable of providing reasonably 
accurate and detailed solutions regardless of the 
complexity of the problem. While the use of 
electrical models to obtain potential and stream- 
line distributions in ordinary flow systems is not 
new, the application of the electrical method to 
problems of the type discussed here does not 
appear to have been previously recognized. 


THE PROBLEM 


The type of problem with which we are con- 
cerned is characterized by the presence of at 
least one boundary of the flow system which is a 
“free surface’ and represents a limiting or 
boundary streamline over which the pressure is 


1M. Muskat, ‘‘The Seepage of Water Through Dams 
with Vertical Faces,” this issue. 


uniform and hence the potential varies linearly 
with its vertical coordinate. Initially the shape 
of this free surface is unknown and its inter- 
section with the outflow boundary surface of 
the system is likewise determined only by the 
complete solution of the problem. Consideration 
of a simple two-dimensional system such as 
shown in Fig. 2 will indicate more clearly the 
nature of the problem. AFDE represents the 
cross section of a permeable dam on one side of 
which water is maintained at height h, and on 
the other at height h,,. It is desired to determine 
the potential and streamline distribution within 
the system with no a priori information regard- 
ing the height of the outflow discontinuity or 
seepage surface BC or the shape of the free 
surface CD constituting the upper limiting 
streamline. 


THE ELECTRICAL CONDUCTION MODEL 


By considering now the construction of an 
electrical conduction model of the system in 
question, it is evident that the inflow boundary 
ED and the outflow surface AB represent 
equipotential surfaces corresponding to the fluid 
heights h, and h,, respectively. Furthermore the 
vertical projection of the terminae of all equi- 
potentials on the boundary BCD must be pro- 
portional to their potentials, the total potential 
difference from D to B being proportional to 
(he—hw). 

Fig. 1 shows schematically the electrical equiv- 
alent of the two-dimensional gravity flow system 
already described. Here a rectangular sheet of 
homogeneous high resistance material A FDE is 
provided with a highly conducting terminal DE, 
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Potential divider 


Fic. 1. The electrical circuit for the electrical model of a 
gravity flow system. 


serving to maintain this boundary at a uniform 
potential e;. Likewise the outflow boundary is 
provided with a similar terminal AB at potential 
és and of height equal to the height of fluid 
maintained at the outflow boundary. The re- 
mainder of the outflow boundary BF is provided 
with a terminal comprising a resistance strip, 
the ends of which are maintained at potentials 
é, and é2 as shown in the figure and having a 
linear potential distribution along its length BF. 
Since the maintenance of a uniform potential 
distribution in this resistance strip, in spite of 
the lateral influx of current from the conductor 
AFDE, would normally require a nonlinear and 
unpredictable resistance distribution, the sheet 
conductor A FDE is made of very high resistance 
material as compared with the terminal strip BF. 
In this manner BF may comprise a resistance 
strip of uniform cross section giving a linear 
potential distribution on passage of current 
through it and yet the perturbation in this 
distribution due to the current feeding into it 
laterally from the conducting sheet will be 
negligible. In the models used by the authors 
the boundary terminal BF comprised a strip 
having a resistance of 0.0123 ohm per cm while 
the resistance of the sheet material was of the 
order of 53.9 ohms per cm (for 1 cm width). 
This contrast is sufficient to eliminate any 
detectable departure from linearity in the po- 
tential distribution along BF. 

The circuit diagram is shown in Fig. 1 and 
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Fic. 2. Potential and streamline distributions in the region 
of flow of a dam with vertical faces. 


requires no discussion except to point out that 
the rheostats R; and Re must be adjusted, with 
the aid of the galvanometer probe circuit and 
with the potential-divider ‘slider set at e;, such 
that the potential at F is exactly the same as the 
potential of the terminal DE, that is, it must 
have the value e,. This equality should be 
checked periodically during the subsequent work 
to make certain that no fluctuations occur either 
through accidental variations in circuit resistance 
or because of the increase in resistance of the 
sheet AFDE incident to the cutting operations 
to be described. A suitable switch or key K is 
provided so that the battery circuit is closed 
only during the short interval required in ad- 
justing the position of the probe point for a null 
reading of the galvanometer. This minimizes 
troubles caused by heating. 

After having set up the electrical system as in 
Fig. 1 it is clear that the known boundary con- 
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ditions of the fluid system have been duplicated. 
Thus is satisfied the requirement that the 
boundaries DE and AB are each equipotentials, 
the difference in values of which corresponds to 
the potential difference between D and B in the 
fluid system. Furthermore, the potential at F is 
equal to that at D and the potential distribution 
along the boundary BCF is linear and propor- 
tional to the height. Evidently, however, the 
shape of the boundary CD, analogous to the free 
liquid surface or upper limiting stream line in the 
fluid system, and its intersection C at the out- 
flow face BF is as yet unknown. In fact the de- 
termination of this boundary is an important 
phase of the final solution of the entire problem. 
In principle, the procedure to be followed in 
determining this boundary in the electrical model 
is actually to cut away the conducting material 
within the area CFD and so to shape the bound- 
ary CD that the potential of any point along 
the boundary CD is proportional to its vertical 
height above point B. Thus by trial is obtained 
the simultaneous solution for the free surface 
and the potential distribution throughout the 
system. Fig. 2 shows clearly that in the final 
result all boundary conditions including those 
for the free surface CD and the “surface of dis- 
continuity”’ BC are properly fulfilled. 

It will now be apparent that proper and easy 
manipulation of the model dictates certain 
physical details regarding its construction. Most 
important is the requirement that the sheet- 
conductor material may be easily shaped by trial 
without otherwise disturbing the system. We 
have found that ordinary Bristol board or other 
heavy firm paper may be given a coating of 
graphite which has the desired moderately high 
resistance, and easily permits the shaping of 
the boundary CD by gradually cutting away 
with a sharply pointed knife or other instrument. 
The use of a moderately sharp needle point in 
firm contact but not penetrating the graphite 
surface serves as a satisfactory potential probe. 
While such graphite material may be available 
commercially we have been able to make such 
sheets by spraying onto the Bristol board, twelve 
to twenty coats of a commercially available 
graphite colloid.? The colloidal solution is made 


= “Aquadag” (concentrated) manufactured by Acheson 
Oildag Company, Port Huron, Mich. 


397 


quite dilute by thoroughly dispersing the colloid 
paste in a suitable quantity of distilled water 
and using it in an ordinary spray gun giving a 
thoroughly atomized spray. Obviously a very 
uniform coating of graphite is required and great 
care is necessary in properly distributing the 
spray. By applying very thin coats and building 
up the final thickness by many individual 
applications, allowing each to dry before further 
application, a very homogeneous coating is 
finally obtained. Before using it in a model the 
homogeneity of the coating should first be tested 
by applying linear electrodes at each end of the 
rectangular sheet and tracing 0.1, 0.25, 0.5, 0.75 
and 0.9 percent equipotentials which should be 
straight lines parallel to the electrodes and pro- 
portionately distributed along the length of the 
conductor. 

The most satisfactory simple method of 
making contact with the graphite sheet is to use 
terminals having very narrow rounded edges 
clamped firmly to the sheet at numerous points 
along the length by suitable clamping blocks. 
The resistance strip used to obtain the linear 
potential distribution (B—F in Fig. 1) is likewise 
a strip resistance material clamped edgewise 
against the graphite sheet. To obtain the 
necessary rigidity the strip is laid edgewise in a 
close fitting slot milled in a strip of Bakelite. 
The strip should project about 1/32 inch along 
the slot and the Bakelite strip should be cham- 
fered so that when clamped in place the potential 
probe may be used right up to the terminal and 
in fact so that the potential distribution along 
the strip itself may be determined. It is advisable 
to determine the uniformity of contact before 
proceeding with the measurements by tracing an 
equipotential very close to each boundary elec- 
trode. Nonuniformity of contact is easily de- 
tected by local distortions of the equipotential 
lines and corrective measures taken by the 
addition or adjustment of clamps. 


PROCEDURE IN OBTAINING FINAL SOLUTION 


In a simple electrical conduction model where 
the boundary conditions are completely defined 
by the statement of the problem and the model 
constructed accordingly, it is a simple matter 
to trace 0.1, 0.2-0.9 equipotential lines or to 
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decrease the interval between equipotentials if 
greater detail is desired. Streamlines may then 
be drawn as orthogonal curves by inspection or 
more accurately by various mechanical aids. Or, 
in the simple models, if a more accurate de- 
termination of streamline distribution is required, 
the problem may be inverted so that the observed 
equipotentials are the streamlines of the real 
system. This inversion only requires that the 
known boundary streamlines be made equi- 
potentials by applying the terminal electrodes 
along such limiting streamlines and the normally 
equipotential boundaries be made insulated or 
“free’’ edges. Clearly the system of equipotentials 
obtained in either case will have the required 
orthogonal relation to the other and an accurate 
mapping of both the potential and streamline 
distribution may be obtained. 

In the present case involving flow under the 
action of gravity, it has already been pointed out 
that the ‘free surface’’ boundary BCD (Figs. 1 
and 2) is unknown. However, the conditions to 
be satisfied along this boundary—that the po- 
tential at any point is proportional to its vertical 
height above B—is defined in the statement of 
the problem. Thus the only additional complica- 
tion introduced in solving the more complex 
problem is the procedure of cutting away the 
conducting sheet until, by trial, the free-surface 
potential distribution requirements have been 
satisfied. Obviously this should be done in very 
gradual stages for once cut below the proper 
point it is impossible to make corrections. It will 
be found that with a little experience combined 
with good judgment obtained thereby, it is 
possible rapidly to approach the final solution. 
It is quite impossible, however, to outline in 
detail more specific directions for this procedure. 
Unfortunately the gravity-flow model is not 
readily ‘“‘inverted’”’ as described above so that 
the streamlines must be drawn as orthogonal 
curves rather than determined by direct ob- 
servation. 


TYPICAL SOLUTIONS 


Since the numerical results are available for 
certain cases of rather simple geometry in the 
analytical paper! already cited, we have selected 
for reproduction in the model, Case I, of that 
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TABLE I. Comparison of analytical and solutions from model. 








ANALYTICAL 
SOLUTION 


EXPERIMENTAL 


DIMENSION SOLUTION 





_— h./c 0.3223 0.322 
panels roe os 
L/e 0.1619 0.157 

c 0.1225 0.123 


h. 








paper. Thus in Fig. 2 which is a direct copy of 
the model itself, dimensions h,, h, and | corre- 
spond, respectively, to dimensions h,/c, h,/c, 
and L/c of Muskat’s example and constitute the 
given boundary dimensions. The final solution 
obtained from the model is shown by the free 
surface CD and the potential distribution ob- 
tained in the manner already outlined. The 
streamlines shown are drawn in as orthogonal 
curves and are only approximate. 

In order to permit a direct comparison between 
the experimental and analytical results, Fig. 3 
was prepared from the data of Fig. 2, the plots 
being directly comparable with those of Fig. 6 
in M. Muskat’s paper.' It should be mentioned 
that the curves showing velocity distributions 
(Fig. 3) are necessarily only approximate since 
they were obtained by measuring the potential 
gradient (reciprocal of the distance between 
equipotentials) at the points in question from 
Fig. 2 and it is clear that the true gradient at a 
given point is not obtainable in this simple 
manner. However, by using the reciprocal of the 
distance between equipotentials at the points in 
question as representing the potential gradient 
and after having established the proper numerical 
coefficient to make the velocity at the inflow 
point E=1.36 (the value obtained analytically), 
the other points on the several velocity dis- 
tribution curves of Fig. 3 were calculated by 
using this same coefficient. 

For convenience Table I shows a comparison 
of the results obtained analytically and those 
from the model where again Muskat’s notation 
is used. While less convincing than a detailed 
comparison of the curves of Fig. 3 and Fig. 6 of 
the reference, h,/c which gives the height of the 
outflow surface is a critical point of the solution 
and the agreement with the calculated value is 
excellent. 

Thus it is evident that the electrical conduction 
models of the type described are capable of 
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Fic. 3. Velocity and pressure distributions for the dam of Fig. 2. 


solving two dimensional problems involving the 
flow of liquids through porous media under the 
action of gravity with an accuracy quite sufficient 
for practical purposes. 

In addition to the information disclosed here 
it is also possible to determine the total flux 
through the liquid system by merely measuring 
the total current through meter J, Fig. 1, with the 
potential divider circuit open, and the maximum 
potential drop e;—é2 across the electrical model, 
after the solution is complete. From these data 
the effective over-all conductance of the system 
may be obtained and since the specific con- 
ductance of the material used in the model is 
easily measured, the over-all conductance of the 
model may be expressed in terms of the specific 
conductance of the medium. This coefficient 
evidently will give the effective permeability of 
the gravity flow system in terms of the specific 
permeability of the. porous medium, and the 
liquid flux through the system may be calcu- 
lated. 

As an illustration of such use of the electrical 
models some interesting results for the case 
shown in Fig. 2 are given. The effective re- 
sistance of the system was measured first with 
the faces AF and ED as uniform potential 
electrodes analogous to a simple linear flow 
system. Next the resistance was measured under 


conditions corresponding to the case of gravity 
flow (Fig. 1) but before section CFD was cut 
away, this being a system which can be solved 
rigorously by analytical methods. Finally the 
resistance was determined for the case wherein 
the free surface CD was cut away and the model 
had the proper potential distribution (Fig. 2) to 
simulate exactly the gravity flow system. 











Fic. 4. Potential and streamline distributions in a dam with 
vertical faces and zero outflow fluid level. 
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Fic. 6. Potential and streamline distributions in a dam whose faces have inclinations of 45°. 


Effective Resistance (ohms) 
Linear 
System 
27.73 
27.96 


Ratio 
Linear/Gravity 
0.633 
0.642 


Gravity System 
(No Free Surface) 


43.78 
43.56 


These results show that for the system in Fig. 
2 having the proper potentials at the boundaries 
AF and ED but without the free surface CD 
cut away, the over-all conductivity is 63.7 
percent of the conductivity of a linear flow sys- 
tem of identical dimensions. Within experi- 
mental error this is in agreement with the theo- 
retical ratio 0.631 obtained by Muskat. 

A comparison of the gravity system with and 
without the free surface has also been made for 
the above case. The results obtained with two 
different models are: 

Effective Resistance (ohms) 
Gravity System 


Without 
Free Surface 


44.1 
43.5 


With 
Free Surface 
44.2 
43.5 


It is apparent that the cutting away of the 
section CFD has a negligible effect upon the 
over-all conductivity and that within the 
experimental error of about 0.25 percent, the 
total flux through the true gravity flow system 
may be assumed to be the same as in the special 
case which may be treated analytically. 

Fig. 4 shows the solution for the case of a 
rectangular dam of permeable material where 
the water level on one side is at D while at the 
outflow face the water level is at A or zero 
height. 

Figs. 5 and 6 are for the case of a permeable 
structure of trapezoidal cross section AFDE 
supporting a water level at D on the inflow face 
but with a zero water level (A) at the down- 
stream side AF. The free surface CD and the 
surface of seepage or discontinuity at the outflow 
face AD and the potential distribution at the 
base AE are of particular interest. 

We wish to emphasize that the particular 
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examples given here are not intended to illustrate 
practical problems but rather were selected 
merely to show the possibilities of the method. 
Evidently a permeable or earth-fill dam would 
not be constructed upon an impermeable base as 
is the case in the present example. For such 
conditions concrete dams are feasible. Further- 
more an earthen dam constructed on permeable 
foundation strata would be provided with an 
impermeable central core as, for example, a 
puddled clay wall extending some distance into 
the permeable substratum. In addition the toe 
or heel of the dam or both might be provided 
with impermeable sheet piling or wall driven 
into the substratum or with aprons extending 
outward from the base of the dam. It is evident, 
however, that such specific conditions may be 
duplicated in the model even to the extent of 
introducing variations in permeability by suit- 
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able preparation of the conducting sheet. Thus 
the conductivity of any portion of the model 
may be made to simulate varying permeabilities 
in the actual case by using suitable masks and 
spraying additional coats of graphite upon those 
portions representing more permeable sections. 
Almost any degree of complexity could be 
reproduced in such models and determinations 
made of the surface of seepage, pressure and 
streamline distribution and important details 
derived therefrom. Sand models have been used 
extensively for such purposes but the electrical 
conduction models are more easily constructed, 
permit higher accuracy, and greater detail in 
the measurements. 

The authors are indebted to Dr. M. Muskat 
for important contributions in this work and to. 
Dr. P. D. Foote for permission to publish this 
paper. 
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The method of hodographs, recently developed by Hamel 
for the analysis of ground-water flow systems containing 
free surfaces, is applied numerically to six cases of the 
seepage of water through dams in which the ratio of the 
inflow fluid head to the width of the dam at the base varied 
from 2.04 to 1.08. In one case the velocity distribution was 

_ computed both along the inflow and outflow faces, as well 
as along the base of the dam. The pressure and velocity 
distributions along the bases and the velocity distribution 
along the inflow faces were found for 5 other cases, in 4 of 
» which the outflow fluid head was zero. The heights of the 
surfaces of seepage and the total fluxes were calculated for 
all six cases. These calculations show that the Dupuit- 


Forchheimer theory of the free surface and velocity distri- 
butions in gravity flow systems is definitely wrong, and 
that the extension of the empirical results of Wyckoff, 
Botset and. Muskat, obtained with radial gravity flow ” 
systems, to linear systems is unjustified. The fluxes, how- 
ever, given by this extension or the Dupuit-Forchheimer 
theory turn out to be in close agreement with those calcu- 
lated by the exact theory. The evaluation of the funda- 
mental integral over the argument of the modular elliptic 
function involved in the calculations is extended beyond 
the range given by Hamel and Gunther so as to include all 
real values of the modular elliptic function. 





INTRODUCTION 


NE of the major outstanding problems yet 
requiring explicit solution in the subject 

of the flow of homogeneous fluids through porous 
media is that of the flow of liquids through such 
media under the action of gravity. This type of 
problem is characterized by the fact that at least 


one segment of the boundary of the region of 
flow is a “free surface,”’ the shape of which is not 
known a priori. All that is known about it is 
that it is a streamline surface over which the 
pressure is uniform. However, in contrast to the 
analogous problem of discontinuous motion—as 
in jets—in the ordinary hydrodynamics where 
the uniformity of the pressure implies a uni- 
formity in the velocity, the only immediate 
consequence following here from the constancy 
of the pressure is that the potential along the 
free surface varies linearly with its vertical 
height above a zero plane. Although the analyt- 
ical problem is that of solving Laplace’s equa- 
tion for the velocity potential ®: 


vV*b=0; b=(k/u)(p+re2y), (1) 


k being the permeability! of the porous medium, 
u the viscosity of the fluid (an effectively incom- 
pressible liquid), y its density, p the fluid pres- 


‘The permeability is defined as the volume of a fluid 
of unit viscosity passing through a unit cross section of the 
porous medium in unit time under the action of a unit 
pressure gradient. For the system of units: cc, cm, sec., 
atmos., centipoise, the unit has been named the “darcy.” 


sure, g the acceleration of gravity, and y the 
vertical coordinate, the difficulties of finding 
solutions satisfying all the boundary conditions 
are still very formidable. 

The first attempt to treat such problems 
analytically was that of Dupuit? who assumed, 
in the case of the radial gravity flow into a well, 
that the liquid moved in towards the wells in 
vertical shells with a velocity proportional to the 
slope of the free surface, the square of which be- 
ing neglected in comparison to its first power. 
Other approximate theories, as those of Bose* and 
Kozeny,‘ did treat the gravity flow problems 
from the point of view of Laplace’s equation, 
but no attempt was made in these to find the 
free surfaces. Rather they were chosen as 
streamlines, found in the solution, of ‘‘reason- 
able” shape, without imposing on them the 
requirement of the uniformity of the pressure. 
The first rigorous treatment was that of Hopf 
and Trefftz> who solved the problem of the 
drainage of a sloping terrain by a ditch dug at the 
top of the sand. Here, however, an indirect 
method of solution was used, the shape of the 
ditch being determined only after assuming a 
simple form for the arbitrary function in the 

2 J. Dupuit, Etudes theorique et pratiques sur le mouvement 


des eaux (1863). 

3N. K. Bose, Punjab Irrigation Research Institute 2, 
No. 1 (1929). 
M eed Wasserkraft u. Wasserwirtschaft 22, 120 
1 . 

5L. Hopf and E. Trefftz, Zeits. f. angew. Math. u. 
Mech. 1, 290 (1921). 
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generalized Schwarz-Christoffel integral for the 
conjugate function transformation. 

Among the experimental studies of the problem 
of gravity flow through porous media, the most 
recent is that of Wyckoff, Botset and Muskat.*® 
The essential conclusion of these authors, work- 
ing with radial sand models, was that while the 
Dupuit formula for the free surface in gravity 
flow systems was incorrect, it was formally 
accurate if the heights of the free surface in that 
formula were interpreted as the fluid heads at the 
base (horizontal) of the flow system. Further- 
more, with this interpretation, the fluxes through 
such systems were also correctly given by the 
equation following from the Dupuit assumptions, 
namely : 


Q= wkyg(h?— Nw) / log re[Tws (2) 


where h., 4» were the fluid heads at the external 
and internal circular boundaries of radii 7., rw. 

These results were also found to apply to com- 
posite flow systems in which a nongravity pres- 
sure differential was superposed on that due to 
gravity, this added differential, however, not 
being great enough to force the fluid to fill out 
completely the sand and thus eliminate the 
gravity component of the flow; that is, when the 
fluid height at 7, exceeded the sand thickness, 
while that at r,, was less than the sand thickness, 
the resultant flux through the system and pres- 
sure distribution at the base were shown em- 
pirically to be given simply by the sum of the 
gravity and nongravity components. These re- 
sults were then generalized—by postulation—so 
as to give Laplace’s equation: 


v*y?=0, (3) 


for the fluid head distribution’ at the base (hori- 
zontal) of every gravity flow system, and a 
general procedure was given for obtaining the 
flow rate through any such system, and also 
for those having both gravity and nongravity 
components. 

The greatest advance, however, yet made in 


' the study of gravity flow systems is that recently 


| ®&R. D. Wyckoff, H. G. Botset and M. Muskat, Physics 
| 3, 90 (1932). 

7 This equation had also been proposed and used pre- 
viously by Forchheimer, with, however, the interpretation 
of y as the height of the free surface above the horizontal 
impermeable base. 
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made by Hamel.* Here, by the use of the method 
of hodographs in conjunction with the theory of 
conjugate function transformations, including 
those defined by the modular elliptic functions, 
a rigorous treatment was given of two-dimen- 
sional systems in which there are not only ele- 
ments of free surface but also of “surfaces of 
seepage,” as will be described below. It is in 
this latter respect that Hamel’s theory makes a 
particular advance over that of Hopf and Trefftz, 
where this additional type of surface element was 
disregarded, in addition to its provision for tak- 
ing into account directly the shape of the outflow 
surface.® 

Unfortunately, however, the numerical calcu- 
lations required to get specific results from 
Hamel’s theory are quite tedius and lengthy. 
The only calculations thus far carried out are 
those of Hamel and Gunther"? in which the shape 
of the free surface for a single case of a problem 
corresponding to the seepage through a dam with 
vertical faces was computed. No results are 
available for the total flow through the system, 
which is probably of greatest physical interest. 
In this paper the results are given of the computa- 
tion of this quantity not only for the particular 
case of Hamel and Gunther, but also for one 
other in which the fluid head at the outflow face 
is nonvanishing, as well as for 4 cases with 
vanishing outflow fluid head. In addition, the 
velocity and pressure distributions have been 
found along the inflow face and base for each 
case, and also those at the outflow face for that 
of Hamel and Gunther. These are compared with 
those which would be implied by the generaliza- 
tion of the results of Wyckoff, Botset and Muskat 
already mentioned. Finally, some additional 
analytical and numerical details are given which 
will considerably lessen the work of studying 
still other cases. 

In order to make clear the principle of the 
computations, a brief outline of Hamel’s theory 


8G. Hamel, Zeits. f. angew. Math. u. Mech. 14, 129 
(1934). 

® Classes of complex variable functions corresponding to 
free surfaces had already been found by Kozeny (Wasser- 
kraft u. Wasserwirtschaft 26, 28 (1931)), but no procedure 
was given for finding the free surface corresponding to 
any preassigned physical system. 

10 The author is very much indebted to Dr. Hamel for 
sending him before publication the manuscript giving the 
detailed calculations for this case. 
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‘will be given first, as applied to the specific 
. problem of the seepage through a dam with 
vertical faces." 


R&ésuME oF HAMEL’s THEORY 


Hamel’s method is essentially that of first 
mapping the original geometrical outline of the 
flow system—the z plane—onto a ‘“‘hodograph 
plane,”’ in which the independent variables are 
the velocity components, and then mapping this 
hodograph plane onto an upper half of a new 
plane, the X plane, the boundary of the flow 
system in the former being transformed into 
the real axis of the latter. The construction of 
the figure in the hodograph plane is obtained by 
the following transformations of the various 
types of boundary segments. 

(a) Impermeable boundary of inclination a to the 
x axis. In the notation of Eq. (1), the velocity 
components (u, v) will be given by: 


OP ad v 
u=——cosa; v=——sina; -—=tana, (4) 
Os Os u 


where ds is a boundary element. Hence for a 
rectilinear impermeable boundary, the hodo- 
graph is a parallel line passing through the origin. 

(b) Constant potential surface—stationary body 
of liquid. Here by definition #=constant, and 
the velocity components at the surface are: 


OP — OP v 
u=——sina; v=+—cosa; —=-—ctna, (5) 
on - On u 


where » is the normal to the surface, so that a 
rectilinear constant potential boundary gives a 
line through the origin in the (u, v) plane normal 
to the real boundary. 

(c) Free surfaces. As by definition the pressure 
is constant along the free surface, © is given by: 


o—cy= C; c=kyg/n, (6) 


where the vertical coordinate, directed upwards, 
has been taken as y. After noting that the free 
surface must also be a streamline of the system, 
it readily follows that along it: 


v+u?+cv=0, (7) 


‘It is assumed throughout that the dam has a length 
that is large compared to its breadth, so that the problem 
is essentially two dimensional. 
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Fic. 1. The z plane representation of a dam with vertical 
faces. 


which is a circle passing through the origin, of 
radius c/2 and center at (0, —c/2). It may be 
noted that c is the vertical free fall velocity due 
to gravity. 

(d) Surface of seepage. At a ‘‘surface of seep- 
age,””"? through which a liquid passes and enters a 
region free of both porous medium and liquid, 
the pressure is constant, although the surface is 
not a streamline surface. From this it may be 
shown that: 


c sin a+u cos a+v sin a=0, (8) 


so that if the surface is plane, its hodograph is a 
straight line normal to the surface and passing 
through the point (0, —c) which is the lowest 
point of the circle of Eq. (7) representing the 
free surface. 

By applying these results to the problem of 
the seepage through a dam, indicated diagram- 
matically in Fig. 1, it is readily found that its 
representation in the hodograph plane is as shown 
in Fig. 2. It is to be noted that whereas in the z 
plane the exact form of the boundary segment 
DC, and hence also the terminus of BC, is in 
principle unknown until the problem as a whole 
as been solved, the boundaries of the system in 
the (u, v) plane are all known except for the 
numerical values of the abscissae of EF’, A’. 

The next step is that of finding another func- 
tion of (u, v) whose actual values on the boundary 


'2 The literal translation of the German “hangquelle”’ 
hardly seems as expressive as the term “‘surface of seepage. 
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Fic. 2. The hodograph plane representation of Fig. 1. 


segments can be preassigned. Such a function 
can be obtained as follows: Introducing first the 
complex potential, 


f(z) =+1¥, (9) 


where ® is the velocity potential and W is the 
stream function, one proceeds to the inter- 
mediate potentials 7, 6 defined by: 


r+i0= —log {—f'"(z)}, 


where since: 


(10) 


(—d(f'(2))/dt_ wu—i be~ix 
dz/dt 





—f"(2)= , (11) 


ut+iv ne‘ 


t=logn/b; 0=w+ x, (12) 


where 7, b are the absolute values of the velocity 
and acceleration and x, w are the inclinations 
of their directions. On the basis of these defini- 
tions it may be shown that the distribution of @ 
on the boundaries is that given in Fig. 2, where 
the angle 8 relating to the free surface has been 
defined by the relation: — 


u=csin B/2cos 8/2; v=—csin?B/2. (13) 


Before mapping the hodograph in the (uw, 2) 
plane on the infinite half-plane it is convenient 
to transform it first into the standard form of 


Fig. 3 obtained by applying the transformation 
function 


q=1+(i/c)(u—w), (14) 


giving the representation of Fig. 3. 
The region B’’ C’’ D” E”’ A” B”’ in the q 
can now be mapped onto the upper half 
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Fic. 3. The q plane diagram corresponding to Fig. 2. 





X plane 
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Fic. 4. The \ plane map of Fig. 3. 


of the \ plane by means of the modular elliptic 
function (gq), this function being so chosen 
that: 


(0) =0; A(1)=1; AGo)=—-; 


M1+io)=+, (15) 


giving the representation in Fig. 4. 
As @ is known on the real axis of the \ plane, 
the complex function @—ir may be expressed 


in the whole upper \ plane by the generalized 
Poisson integral: 


1 pt 0(t)(At+1)dt 


6—ir = —iro+— 


ri Jo (t—d)(1+2) 





where 79 is an arbitrary constant. Denoting now 
the values of \ corresponding to E’’’ and A’’’ by 
b and a, evaluating the integral of Eq. (16), 
and combining with Eq. (10), one finds: 


(\—b)(a—d)AT} 
— 


3 7 Bd) 
xX =e ——d e 17 
exp +— f t—xr | (17) 





s"@)=-| 
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As \=X(q) is a known function of u—iv= —f’(z), 


it follows that f’’(z) = F(f’(z)), and: 


d(u —iv) 
z= -f ——= f er*td(u—iv) 
F(u—iv) 
h-1 , 
-cf| | 
\(a—d)(A—}) 
3" BCL) 
exp| -— f — at |d(u—io), (18) 
To t—xX 


where C=e, and the sign of the radical for 
real \ is to be taken as that of e~‘®. Thus, in 
principle, the problem is solved, as by Eq. (18), 
z is given as the function of u—iv, and hence the 
velocity components as functions of z. 





SPECIFIC CALCULATIONS 


The physical content of the analytical solu- 
tion contained in Eq. (18) may be obtained by 
the evaluation of the geometrical dimensions of 
the flow system as implied by the constants a 
and b, the determination of the shape of the 
free surface, the computation of the velocity 
distributions along the boundaries, the calcula- 
tion of the total flux through the system, and the 
mapping of the potential and streamline distri- 
butions in the system. The specific expressions 
giving these results are readily seen to be the 
following equations, taking note of the values 
of @ indicated in Fig. 2, and setting the constant 
C equal to unity: 

(a) Velocity distribution (x=x(u)) along EA, 
its length LZ, and the fluid head distribution 
head h(x): 


u(A) 
edu, 
u(E) 


x(EA)= edu; L= 


u( BE) 


(19) 


hx) =he— f udx=he— f uedu. (20) 
0 (BE) 


u 


(b) Velocity distribution (y.=yw(u)) along 
the outlet face AB, and the outlet fluid height 
hw: 


fos) 


u 
ve f edu; he= f 
u(A) u(A) 


edu. 


(21) 
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(c) Velocity distribution (y,=y,(u)) along the 
surface of seepage BC, and the length h, of BC: 


u 0 
y= -{ edu; h,= -{ evdu. (22) 
u( B) 


@ 


(c) The shape of and velocity distribution 
(x=x(u, v)), (yv=y(u, v)) along the free surface 
CD: 


2—2.=x—L+1(y—hw—hs) 


. c 78 
= f esera(u— in) =f e7— ‘8/28, (23) 
(e) Velocity distribution (y.=y.(u)) along 
the inflow surface, and the inflow fluid height h,: 


u 0 
= -f edu; h.= -f edu. 
(B) (EB) 


(24) 


u u 


A convenient check between the computation of 
h. and h, is obtained from Eq. (20) in the 
relation: 


(25) 


L u(A) 
he—he= f udx= f ue"du. 
0 u(E) 


(f) Total flux, Q, through the system: 


0 


he hw 
o=f udy, = -{ uerdu= UudYw 
0 (EB) 0 


u 


he co 0 
+f udy,= { uerdu— f uetdu. (26) 
0 (A) 


u @ 


(g) To find analytically the potential and 
streamline distributions in the interior of the 
system would require the general evaluation 
of the integral of Eq. (18) so as to first get the 
internal velocity distribution and then (x, y) 
by integration. However, the most rapid and 
convenient procedure would probably be to find 
(x, y) by standard numerical or graphical 
methods for solving Laplace’s equation after 
having determined the shape of the free surface | 
by Eq. (23) and the values of ® along EA by, 
Eq. (20).% 

134 still more rapid empirical method of finding the 
shape of the free surface and the internal potential dis- 
tribution, and one applicable to systems’ for which the 
analysis would be quite intractable, is that based on the 


electrical analogy of the fluid flow through a porous medium 
(cf. R. D. Wyckoff and D. W. Reed, this issue). 
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On the other hand, from a practical point of 
view, the most important feature of the solution 
is the relation of the total flux, Q, through the 
dam to the inlet and outflow heads h,., hy. We 
shall, therefore, stress mainly this phase of the 
solution, and present the other results as they 
are derived incidentally to the calculation of Q. 

In order, however, to get any numerical re- 
sults at all, it is necessary to determine the value 
of the integrand in Eq. (18). Here, it must first 
be noticed that whereas the modular elliptic 
function A(g) has been defined by Eq. (15), its 
properties are usually given in terms of the 
function &? = @(q) defined by the conditions: 


Q(0)=1; O(1)==; 


O(ix)=0; O(1t+ix)=0, (27) 


from which it is seen that 


A(q) =1—1/0(q)=O(1—-1/9), (28) 


which can be expanded by the rules already 
established for the function 9. This relation 
affords the link between the variable \, entering 
in the integrand e’ of Eq. (18) and the variable 
of integration u—iv=—ci(q—1), and thus also 
permits changing the variable of integration 
from u—1v to X. 
Finally, one must evaluate the integral: 





—2r 1 B(t)dt 


use being made of the relation: 


407 


3 "Bit 
an)=—— f , 


(29) 
2r t—X 


This may be done by noting that the variable ¢ 
in Eq. (16) really represents the real axis of the 
\ plane on which the gq plane of Fig. 3 is mapped 
by the function A(g). Hence: 


t=X(q) =d(3+e'8/2) = O(1 —1/¢) 


=O(i tan B/2)=k?, (30) 
from which 8 can be computed as a function of ¢. 

The evaluation of a(d) for 1.05=A=20 and 
—20=A= —0.01 has been carried out by Gun- 
ther, using direct graphical integration. How- 
ever, in evaluating the height h, by Eq. (24) and 
the flux Q by the first of Eq. (26)—the simplest 
methods if one is to avoid computing in detail 
the shape of the free surface—it is necessary to 
know the values of a for \ much closer to 1, as in 
the variable of integration u, \=1.05 corre- 
sponds to u/c=0.5423 whereas the integration 
has to be extended to u=0 (cf. Eq. (24)). 
Since \=1 is a singularity of a(A),—a(A~1) 
~3 log log (A—1)—graphical methods are in- 
accurate. An analytic expansion for very small 
n=A—1, was, therefore, derived as follows: 
From the definition of @ it is seen that: 


, One 0-1 B(#) 
~y wee Eee 
0. 0 t+n 0.01 


B(t)=x—B(1—2), 


(32) 


and the expansion of (¢+ )~! in the first integral being replaced by ¢—. A graphical integration of this 
integral gave the value 8.926. In the remaining integral, 8(¢) was expanded as: 


TT T 


3 





B(t) 1 16 t\\ = 
—=tan"! -( tog —--) =-— 
2 T t 2 


+ oe 
2 log 16—t/2—logt 3(log 16—#/2—log #)* 


(33) 


the argument of the arctangent being accurate to terms in #. By replacing log 16—¢/2 by: log 16 
—0.0025 =log 6, the integral in question reduced to a series involving the integral: 


0.01/65 
1=-f 
0 


and its derivatives with respect to log 6. 


4 In the manuscript already cited, reference 10. 


dz 





(2+n/6) logs 


(34) 
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I was then evaluated by dividing into the ranges 0/6 and 7/5—0.01/6, and expanding the fac- 
tors (s+ 7/6)~'. The final result was found to be: 


—2r 0.01 7 7° 





1 
+0(100n) +0( —— 
log® 9/6 


which is accurate to about 10~* up to 7=10~*. 

Fig. 5 gives a plot of e* as a function of 
—log (A—1), the scale and coordinates being 
chosen so as to permit a representation of 
e* to reasonably small values of u, the value 
of —log (A—1)=27.63 corresponding to u/c 
=0.1033, although A\—1=10-". The scale be- 
comes quite unsatisfactory for \>2; but for such 
values, e* is listed by Hamel and Gunther. 

For large values of \, e* is most conveniently 
obtained by expanding (¢—d) in the integrand of 
a and evaluating the resultant integrals graphi- 
cally. In this way it is found that: 


e% = e3/4d+0.3286/ 2+ 0.2029/AF +--+ | (36) 
- Finally, the values of e* for \<0 can be ob- 
tained from those for \>0 by the formula: 

e*(k) =[ —A/(1—A) J¥e*(1—A), (37) 
as may be readily shown by changing the vari- 
able in the integral for a and then applying 
Eq. (32). 

The connection between \, the argument of 
e* and hence of e’, and the variable of integration 
u occurring in Eqs. (19) to (26), may be‘made as 
follows: From the theory of the modular elliptic 
functions ©0(£)=?, it has been shown" that 0 
may be inverted as: 


aK’ _ 
rit = ——— =G(k?) 
K 


k ® 13k 32_ 
=2 log -+—+——_+—k*+---, (38) 
4 2 64 192 


where K, K’ are the complete elliptic integrals 


Cf. H. A. Schwarz, Formeln und Lehrsatze zum Ge- 
brauche der Elliptichen Functionen, §40. 


' B(z) 
a= rn log (40.01) + f —dz+27 log log ——+ _ 
0.01 2 6 3(log 1006)? 10(log 1006) 


n 
— 27 log log— 
6 


) =7.1927-25 log {log 1/n+2.7701}, (35) 





of the first kind, with modulus &™ and k’, re- 
spectively, the individual values as well as their 
ratio being listed in complete form by Hayashi.!” 
Noting now the relation between \ and © as 
given by Eqs. (28) and (30), and the transforma- 
tion formulas between the function O[ (m’ +n’ €)/ 
(m+né) ], m, m’, n, n’ being integers of mod. 2, 
and ©(£),'!* we readily find that for: 


—x<dA=0, 


1—X K 


” 1 aK’ 1 
kP?=O(qg)=——_-;_ riq= -*_-c(—). (39) 
r 


so that, noting that for \<0, v= —c, g=iu/c, 


u/c=K'/K=(-—1/7r)G(1/(1—A)); 


Similarly, for: O0=A=1, k® may be taken as 
1—X or X, so that: 


nig aK’ 
—— = -— =G(1-)); 
1-—g K 


(¢q—1) 1K’ 
si = ——=G(a); 
q K 


(41) 


16 The modulus of the elliptic integral K is here taken 
as k in order to avoid confusion with the symbol & already 
used to denote the permeability of the porous medium 
(ef. Eq. (1)). 

17K. Hayashi, Tafeln der Besselschen, Theta-, Kugel-, 
und anderer Functionen, 1930. 

18 Schwarz, reference 15, §§33, 54. 
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-log (A“D) 


Fic. 5. Plot of e*, defined by Eq. (29), as a function of —log (A—1). 


and, for: 1=\= ~~, k® may be taken as 1/) or 
(A—1)/X, so that: 


u K’ -1 gl 
ole 


cf hUO« r \A 


c K’ -1 sr-1 
~-—=—6(——), 
uK 7 r 


With the aid of these relations the actual 
integrations in Eqs. (19) to (26) may be most 
readily carried out for finite values of \ and for 
values not too close to singular points of e’ by 
plotting e’ as a function of u/c and integrating 
graphically. In evaluating the integrals for 
infinite values of u or i, the integration in u 
beyond 2.8604 c(\=500) was carried out ana- 


(42) 


lytically in the variable \ with the aid of the 
approximations: 


e7=[(A—1)/A(A—) (A—D) Jhe* =1.001/X, 

(43)!9 
u=(c/r){—1/2A+log 16}; du/dd Sc/mh, 
the latter two equations following from the first 
of Eqs. (42). For the integrations at the singular 
points, as \=b, the evaluations were again 
carried out analytically from }\=) to A\=d plus 
0.01 or 0.05 by replacing everything in the 
integrand except (A—b)-! by its average over 
the interval of integration. In spite of the 
graphical and numerical character of the work 


19 These approximations are the same as those used by 
Hamel and Gunther. 
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it is felt that the errors in the final results should 
not exceed 1 percent. 

For the four cases computed with vanishing 
outflow fluid level, which correspond to a=, 
the factor (a—X) in the above expressions was 
omitted. That this is necessary may be readily 
shown by repeating the analysis leading to 
Eq. (18). Eqs. (19) to (26) remain unchanged 
except that u(A) must be replaced by ~ =u(B). 


RESULTS 


The calculations were carried out in the man- 
ner described above for the six cases defined by 
the values of a and 3b, given in the second and 
third rows of Table I. The results are listed in the 
succeeding rows. The values of Q in this table 
have been computed by means of the first of 
Eq. (26). An indication of the accuracy of the 
calculation of this quantity is given by the fact 
that for case I, Q/c? was found to be 0.2988 by 
the 3rd and 0.2979 by the last of Eq. (26). 
The values of h./c or (h.—h,)/c were also 
checked by means of Eq. (25) in all cases except 
II, where Eq. (25) was used to obtain h,,/c, the 
difference between the values given by Eq. (24) 
and Eq. (25) always being less than 1 percent. 
The value of Q/c? are those computed by the 
formula analogous to Eq. (2), applying to 
linear flow, namely: 


O=(kyg/u)(h2—hy?)/2L, (44) 


which follows both from the Dupuit assumptions 
and the generalized interpretation of Eq. (3) as 
given by Wyckoff, Botset and Muskat. 

The shape of the free surface (calculated by 
Hamel and Gunther) and the velocity distribu- 
tions along the inflow and outflow surfaces as 


TABLE I. 








Ill IV 

10 5 

x io) 
0.67204 0.8715 

0 0 


0.4300 0.5192 
0.3293 0.4843 
0.6871 0.7828 


0.6865 0.7841 
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well as along the base are shown for case I in 
Fig. 6. For cases II-VI, neither the shape of the 
free surfaces nor velocity distribution along the 
outlet faces were computed. However, the 
heights at C of the surface of seepage were calcu- 
lated, and at least a rough idea of the actual free 
surface should be given by the free surfaces 
arbitrarily drawn in the Figs. 7, 8 and 9 for 
cases II, III, and VI, following the general 
features indicated by that for case I. The velocity 
distributions along the inflow faces and the 
bases for these other cases are given in Figs. 7 to 
9.2° These figures also contain the pressure distri- 
bution (calculated by Eq. (20)) along the base, 
and for comparison those, p, implied by the 
solution of Eq. (3), namely: 


p?= —(h?2—h,*)x/L+h2. (45) 


DISCUSSION 


Although the above results are essentially 
numerical and include only six cases of a single 
type of gravity flow problem, they nevertheless 
imply some important consequences which 
should be of significance with respect to other 
gravity flow systems. First, however, one may 
note some of the features of the free surface, for 
case I, shown in Fig. 6. Thus, on noting from 
Eq. (13) that the inclination of the free surface is 
simply 6/2, it follows from Eq. (23) that the 
slope of the free surface decreases monatonically 
from zero at the inflow face at D to a negatively 
infinite”? value at the surface of seepage at C. 
Furthermore, Fig. 6 shows that the actual heights 
of the free surface may be very much different 
from those computed by Eq. (3) with y in- 
terpreted as the fluid height—which would be 


20 The same data as indicated in Figs. 7 to 9 were also 
computed for cases IV and V, but their similarity to those 
shown makes their reproduction unnecessary. The dis- 
crepancy between the dimensions shown in Fig. 9 and those 
given in Table I is due to an error in the calculations dis- 
covered after the figure was drawn and sent in. 

*1 The statement by Wyckoff, Botset and Muskat 
(reference 6) that the slope can never exceed 45° was 
based on the erroneous assumption that the horizontal 
velocity component at the free surface would increase 
monatonically in approaching the outflow surface. How- 
ever, this does not invalidate their conclusion that the 
horizontal velocity component at the free surface can 
never exceed half of the vertical free fall velocity, c. It 
may also be noted that the free surface will always become 
tangential to the outflow face as it joins the surface of 
seepage, as long as the inclination of the outflow face 
equals or exceeds 90°. 
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Fic. 6. Free surface, velocity, and pressure distribution for a dam defined by: b=5, a=10, 
case I. a/c and # are the pressure and velocity distributions implied by Eq. (45). The velocities 


along EA have been reduced by a factor of 10. 
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Fic. 7. Velocity and pressure distributions for a dam defined by: b=2, a=5, case II. a/c, 
and # are those implied by Eq. (45). The free surface was drawn in, without calculation, so as to 
terminate at C. The velocity distributions along EA have been reduced by a factor of 10. 


given by the curve marked p—as originally 
suggested by Dupuit and used at length by 
Forchheimer. Especially in such cases as III-VI 
where the outflow fluid head is zero, Eq. (3) 
would also give vanishing fluid heights near the 
outflow surface, whereas the actual outflow 


heights as given by h, may become quite large 
fractions of h, (for III, h,/h.=0.64). Although 
these discrepancies rapidly decrease as one re- 
cedes from the outflow surface, and should also 
be smaller at distant points in radial systems 
than in linear systems, it does seem difficult to 
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u/c x/c 
Fic. 8. Same as Fig. 7, with b=10, a= ~, case III. 


x/c 


Fic. 9. Same as Fig. 7, with b=1.2, a= ~, case VI. 





SEEPAGE OF WATER THROUGH 


maintain that the Dupuit-Forchheimer predic- 
tions of the shapes of free surfaces are more 
reliable than reasonable guesses, except at points 
where the slopes of the free surfaces are definitely 
known to be small. 

These serious discrepancies must evidently 
result from erroneous assumptions. That the 
Dupuit assumption of ‘shell flow’’—uniformity 
with depth of the horizontal velocity—is falla- 
cious, is clearly shown by the velocity distribu- 
tion curves of Figs. 6 to 9. Obviously the velocity 
variation from 0 at D to the maxima at E can 
be put poorly approximated by a uniform veloc- 
ity distribution.” Similarly, the variation of u/c 
along AB in Fig. 6 from a minimum at A to « 
at B is hardly suggestive of any essential uni- 
formity, even though the region of indefinitely 
large growth of the velocity is quite limited. 
Finally, the surface of seepage BC, which is 
entirely neglected in the Dupuit-Forchheimer 
theory, also possesses a range of infinitely high 
velocities. In fact, in case I, the flow through 
BC happens to be almost exactly equal to that 
through AB, and in such cases as III—VI where 
h.=0, the total flow must, of course, pass 
through BC. It is, therefore, clear that the Du- 
puit-Forchheimer basic assumptions and inter- 
pretation of Eq. (3) as a differential equation for 
the free surfaces of gravity flow systems must 
be used with great caution if entirely erroneous 
results are to be avoided. 

The interpretation of Wyckoff, Botset and 
Muskat® of Eq. (3) as a differential equation for 
the pressure distribution at the base (horizontal) 
of a gravity flow system also seems to be subject 
to considerable error. For, as will be seen from 
Figs. 6 to 9, the pressure distributions, p, im- 
plied by this interpretation, differ quite ap- 
preciably from that given by the rigorous theory. 
The discrepancies are particularly marked near 
the outflow faces if h,,=0 (Figs. 8 and 9), where 
Eq. (45) leads to a velocity that becomes 
infinite as x}, x being the distance from the 
outflow end, whereas the exact analysis re- 
quires a logarithmic infinity for u. As a whole, it 


22 Tn fact, if the Dupuit method is followed to its logical 
conclusion, and one uses for the average horizontal velocity 
along the vertical section its actual value at the free 
surface, the absurd value of zero flux would be obtained 
when this procedure is applied to either the inflow or 
outflow surface. 
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will be noted that the actual pressure distribu- 
tions at the base are much more linear than is 
implied by Eq. (45). 

While this result may appear surprising in 
view of the close agreement found by Wyckoff, 
Botset and Muskat between the predictions of 
Eq. (3) and the pressure distributions obtained 
from actual sand model experiments of radial 
flow systems, the reason undoubtedly lies in the 
differences in geometry of the systems. Thus in 
the radial system, where the effective resistance 
is largely concentrated near the outflow surface, 
a distortion of the geometry in this region, such 
as is caused by the shrinkage of the actual height 
of sand through which liquid is flowing due to 
the depression of the free surface,”*> might well 
result in a marked modification of the pressure or 
fluid head distribution (even at the base) from 
that in a strictly nongravity radial system, 
namely: 


(he—hw) log r/tre 


h 





+hy. (46) 
log re/Tw 


And indeed, it was found there that the ob- 
served pressure distribution could be expressed 


by: 


(h22—h,”) log r/r 
h?= 





+h,,’. 


(47) 


log re/Tw 


However, in the case of the linear flow system, 
the region near the outflow surface is no more 
sensitive to distortions in geometry than any 
other part of the sand, and is moreover not at 
all as sensitive as that in a radial flow system. 
Hence the effect on the pressure distribution of 
the depression of the free surface will necessarily 
be smaller, and it will be only slightly modified 
from that in the strictly linear system, except 
where h,,=0 and infinitely high velocities are 
thus necessarily induced. It is, therefore, not 
unreasonable that Eq. (45), which is based upon 
the extrapolation of the effect of gravity ob- 
served for radial flow—the change from Eq. (46) 
to Eq. (47)—should be found to give too large a 
modification in the case of linear flow. 


23 The shrinkage of the cross-sectional area due to the 
depression of the free surface has the same type of effect 
as decreasing the permeability of the sand. 
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Although these considerations seem to leave 
Eq. (3), even with the interpretation of y as a 
fluid head, restricted essentially to radial systems, 
it, strangely enough, retains its significance in 
providing a satisfactory method for computing 
the flux in gravity flow systems, a knowledge of 
which is of even greater practical importance in 
the discussion of gravity flow systems than the 
shape of the free surface. For not only was Eq. 
(2), which is a consequence of Eq. (3) for radial 
flow, verified directly and empirically in the 
course of the establishment of Eq. (47), but a 
comparison of the last two rows in Table I shows 
that the linear flow implication of Eq. (3) also 
checks satisfactorily, for practical purposes, the 
correct values (Q/c?). Thus, while the calcula- 
tions carried out here definitely disprove the 
validity of Eq. (3) for giving either the fluid 
height distributions or the fluid head distribu- 
tions (at the horizontal lower boundaries) in 
general gravity flow systems, they at the same 
time give strong support to the equation when 
used simply as a means of computing the flux 
through such systems, although it is difficult 
to find any rational explanation for the apparent 
cancelation of all the obvious errors in the solu- 
tions for the y’s in such a way as to nevertheless 
give close approximations to the resultant fluxes. 
(See, however, the ‘‘note added in proof.’’) 

The method for finding the flux from the 
solutions of Eq. (3) is simply contained in the 
equation: 


dy kyvyg roy’ 
—ds=— | —ds, 
2u on 


kyg 
Q=—I|y 
be on 


(48) 


where the integral is extended over either the 
external or internal boundary of the system, 
with normal n, the solution y? of Eq. (3) being 
that which assumes the given values of h,,”? and 
h.2, the squares of the fluid heads, at the outflow 
and inflow surfaces. This procedure will give the 
same numerical results as that of Dupuit,™ 


* It is to be observed that derivations of Eqs. (2) and 
(43) by the Dupuit-Forchheimer theory involves the 
assumption, in addition to those mentioned in the intro- 
duction, that there is no surface of seepage. A consistent 
application of that theory would require that h, in these 
equations should be replaced by h»+h., and hence lead to 
results which do not agree either with the empirical or 
rigorously calculated values for the fluxes Q. 


MUSKAT 


although its only justification is that it gives 
values of Q which agree quite closely with those 
found either by direct experiments in the radial 
flow case or exact calculation in the simple case 
of linear flow.” 

Summarizing, then, it appears that the calcu- 
lations carried out here not only show the nature 
of the correct solution of two-dimensional gravity 
flow problems, but moreover they for the first 
time conclusively prove the inadequacy of 
previous approximate treatments and empirical 
generalizations with respect to the shapes of the 
free surfaces or pressure distributions within the 
gravity flow systems. Finally, they show at the 
same time that these earlier formal approxima- 
tions do lead to predictions of the flux through 
gravity flow systems—at least for radial and 
linear systems—which are accurate to within one 
percent. 

Note added in proof: A closer examination, 
made since this paper was submitted, of the 
physical significance of the second method, men- 
tioned in footnote 25, of deriving Eqs. (2) and 
(44), has revealed the real reason for their success 
in giving a very close approximation to the flux 
through gravity flow systems. Thus if we replace 
the linear gravity flow system of Fig. 1 by one 
with exactly the same boundary conditions along 
the faces but without an upper free surface 
boundary, i.e., with the conditions: y=0, h, : 
db/dy=0; x=L, O=y=h.: b=ch.; x=0, 
O0=yS=h,: &=chy, and h,=y=h.: &=cy, 
(c= kyg/p), it is readily found that: 


oa —~__.(h 3— 5 3 
® ap, Met hw?) + ap (he hy?) x 


2che $1 (—1)"—cos(nthw/he) | 


ney. pun(L—x 
mT mm sinh (nxL/h,) 


) 
+ he hs ss 





cos 


25 [t may be noted that both Eqs. (2) and (44) which 
follow from Eq. (48) might be intuitively constructed as 
rough approximations on the basis that they are equivalent 
to the values of Q that would be obtained from non- 
gravity flow systems of fluid head differential (h.—hw) 
and “average thickness” (h.+h,)/2, although the latter 
gross a would become immediately question- 


able when the existence of the surface of seepage, and 
hence outflow height of h~+4., is recognized. A somewhat 
more satisfactory mode of intuitive approximation could 
be based on the observation that the actual average 
potential over an outflow surface, of total height A, 
which is faced by a free liquid reservoir to the height 
hw, is ®Pw=kygihe+h.)/2hmu. Hence assuming a strictly 
nongravity flow through the total sand thickness h,, 
but with an effective potential difference of ®,.—?, 
=kyg(h2—h.?)/2heu, the same Eqs. (2) and (44) will 
again be obtained. 
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The flux corresponding to this ®, or: 
S,'(ab/dx)dy is evidently given by Eq. (44). 
Now as the real gravity flow system differs from 
that represented by the above ® only in the 
corner near x=0, y=h,., which lies above the 
actual free surface, but which is obviously a 
region of very small flux density, it is clear that 
one may justifiably expect that the flux in the 
approximation system (without the free sur- 
face), ie., Eq. (44), should differ but slightly 
from that in the physical gravity flow system. 
And indeed this is verified by reference to the last 
two rows of Table I. Of course, it is evident from 
this rigorous derivation of Eq. (44) from the 
approximation system, that it should give values 
for the flux that are somewhat higher than those 
passing through the gravity flow system with its 
free surface. That the difference is, however, 
extremely small is shown by the fact that within 
the accuracy of the calculations (one percent) 
there appears to be no difference at all between 
the last two rows of Table I (the excess of Q/c? 
by one percent for case I must be due to numeri- 
cal errors in the use of the first of Eqs. (26) for 
this case, as the last two of Eqs. (26) gave an 
average value even less than Q/c?). In fact, this 
is conclusively verified by the direct electrical 
model measurements of Wyckoff and Reed 
(cf. preceding paper) which showed, for case I, 
that the conductivity of the approximation 
system exceeds that of the real gravity flow 
system by no more than 1/4 percent. That the 
Dupuit-Forchheimer theory gives Eq. (44) must, 
therefore, be considered as entirely fortuitous, 
while its success as an approximation formula is 
to be attributed to the fact that it is the correct 
value of the flux through a system which, from 
a priori considerations, could really have been 
expected to give a very close approximation to 
that in the physical flow system. 

Exactly the same type of approximation 
system for the case of radial flow (with circular 
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boundaries: r=frw, 7.) gives for the internal 
potential distribution: 
c(h2—h,?) r 
2h, log (re/Tw) °8 re 
+ 2che$ {(—1)"— cos (nmhw/he)| U(anr) cos (nry/h.) 
7 1 n? U (ant w) : 


where: U(a,r) = Io(nar/h.)Ko(nrr./h.) 
—I(nxr./he)Ko(nxr/h.), 


and a flux identically the same as Eq. (2). Here, 
too, the error involved in retaining in the flow 
system the region above the actual free surface 
should be very slight, as was indeed found to be 
the case in the original sand model radial flow 
experiments (cf. reference 6), the fluxes in which 
were reproduced by Eq. (2) with an accuracy 
well within the experimental errors. 

Finally, it may be noted that the above ap- 
proximation systems not only give very close 
approximations to the flux in the physical 
gravity flow systems, but the pressure distribu- 
tions at the bases have been found to be also 
satisfactorily reproduced by Eqs. (49) and (50). 
This is, of course, to be expected, as the mere 
presence of the free surface should not appre- 
ciably affect the potential distributions in parts 
of the flow system that are distant from it. This 
method of replacing the physical gravity flow 
system by those having the same boundary 
conditions, without the free surfaces, thus 
affords a relatively simple and fairly accurate 
means of studying all the significant features of 
the type of gravity flow system considered here 
excepting those directly related to the free 
surfaces themselves. 

The author is indebted to Mr. R. D. Wyckoff 
for helpful discussions of this paper, to Mr. M. 
Meres for assistance with the calculations, and to 
Dr. P. D. Foote for permission to publish the 
paper. He also again wishes to thank Dr. Hamel 
for permitting him to see the calculations of Mr. 
E. Gunther before their publication. 
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With a new type of especially stable electrolytic capaci- 
tor, the dielectric properties of anodic layers are measured. 
The investigation shows that variations of the dielectric 
properties may be divided into two classes: lagging and 
nonlagging variations. Lagging variations, which include 
variations of the capacity and power factor with the opera- 
ing voltages, are especially characteristic of anodic layers. 
They arise from the fact that the structure (organization) 
of layers adjusts itself with changes of the operating condi- 


tions to a new dynamic equilibrium, within a short transi- 
tion period. On the other hand, in regard to variations of 
the capacity and power factor with frequency, which are 
nonlagging variations, anodic layers are similar to other 
composite dielectrics. An oscillographic study of harmonics 
is made which verifies the lagging nature of the changes of 
capacity with the operating voltages (a.c. and d.c.) and 
demonstrates for the first time the existence of a non- 
lagging variation of the a.c. conductivity with voltage. 





INTRODUCTION 


IGHLY insulating anodic layers formed by 

an electrolytic process on aluminum are the 
dielectrics in the usual electrolytic condensers. 
Although these layers find extensive application 
in practice, the information which appears in the 
literature regarding their dielectric properties has 
many inconsistencies and contradictions. This is 
due partly to the type of experimental measuring 
methods which have been employed in earlier 
investigations, and partly to the lack of uni- 
formity and stability of the specimens used. 
However, work done at this laboratory! has 
shown that anodic layers of a strictly defined and 
stable character can be made. Subsequently, 
by using such layers in conjunction with proper 
electrolytes and adopting new measuring meth- 
ods, it has been possible to make a satisfactory 
study of their dielectric properties. This study 
is the subject matter of the present paper. 

An anodic layer is fundamentally of a dynamic 
nature, and it is therefore important in a sys- 
tematic classification of its dielectric properties 
to differentiate carefully between those effects 
which are instantaneous and those which occur 
only after an appreciable time lag. A lack of this 
distinction is responsible for much of the con- 
fusion which exists in the literature on the a.c. 


properties of anodic layers, and the way in which . 


they are correlated with the d.c. properties. For 
example, if a d.c. potential is applied across an 
anodic layer, a current will flow which gradually 
decreases with time. However, the curve of 

1 J. E. Lilienfeld, L. W. Appleton, W. M. Smith and 


J..K. Nieh, “Studies of Fully Organized Layers on Alu- 
minum,” Trans. Am. Electrochem. Soc. 61, 531 (1932). 


current against time depends not only on the 
applied potential and on the formation voltage of 
the layer, but also on the past history of the 
layer. Thus the conductivity of an anodic layer 
is subject to continual change, depending on 
the conditions under which the layer is operated. 
(Deterioration effects.2) Such changes in con- 
ductivity are examples of what we call lagging 
effects. They occur gradually with time and 
indicate permanent changes in the anodic layer. 
On the other hand, if we measure the capacity 
of an anodic layer at 60 cycles and then at 1000 
cycles, we find a difference in the two capacitance 
values. This, however, does not indicate a change 
in the layer, because if the capacity is immedi- 
ately thereafter measured at 60 cycles, it will be 
substantially the same as before. The variation 
of capacity with frequency is an example of a 
nonlagging effect. 

The properties of dielectrics as usually meas- 
ured are of a nonlagging character. Lagging 
changes are, for the most part, peculiar to 
anodic layers, being due to their dynamic nature. 
Now, in order to obtain a coherent series of 
data on the dielectric properties of anodic layers, 
it is imperative to reduce all lagging changes to 
the utmost attainable minimum. Previous work 
of this laboratory has shown that lagging changes 
become more prominent as the concentration and 
mobility of the active OH~ ions in the electrolyte 
increase. This is true for both the electrolyte of 
formation and the electrolyte of operation of the 
layer. Therefore, the bulk of the experiments 


2 J. E. Lilienfeld, L. W. Appleton and W. M. Smith, 
“The High Voltage Anodic Layer on Aluminum,” Trans. 
Am. Electrochem. Soc. 58, 225 (1930). 
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here presented have been conducted with layers 
formed in electrolytes of low OH™~ concentration 
and operated in extremely viscous electrolytes of 
high specific resistance and low active OH™~ con- 
centration. These latter electrolytes will, in what 
follows, be termed “superviscous electrolytes.”’ 
Their specific resistance is greater than 100,000 
ohms/cm*. In order to keep down the effect of 
this high specific resistance in the power factor 
of the condenser as a whole,’ extremely thin 
coatings of electrolyte have been used, according 
to a technique developed at this laboratory for 
practical condensers, U. S. Patent No. 1,986,779. 
With the use of such specimens, the dielectric 
properties of anodic layers have been investi- 
gated. In addition to the usual capacity and 
power factor measurements, a study has been 
made of the harmonics generated, to yield a clue 
to the law of variation with voltage, where such 
variation exists, of nonlagging changes of the 
capacity, power factor, and a.c. conductivity. 

In experiments where a large percentage of 
harmonics has been desired, anodes formed and 
operated in decidedly alkaline aqueous electro- 
lytes have been used, because they accentuate 
harmonics. As already mentioned, in such anodes 
lagging changes are also accentuated. 

Although this paper is concerned with the 
purely physical investigation of anodic layers, 
the superviscous specimens used in it are also of 
interest because of the application of the prin- 
ciples embodied in their design to commercial 
electrolytic condensers in order to secure a high 
degree of stability and a relatively low power 
loss. This is desirable for a number of purposes, 
e.g., for continuous operation on a.c. circuits or 
on d.c. voltages which vary within very wide 
limits; inasmuch as less stable electrolytic con- 
densers will be subject to a progressive deteriora- 
tion if operated under these conditions. 


ELECTROLYTE POTENTIAL 


An anodic layer can be studied by the usual 
methods designed for the investigation of good 
dielectrics only so long as it does not conduct 
appreciable current in either direction. This 


3 The power factor of an electrolytic condenser is due to 
the resistance of the electrolyte as well as the loss in the 
anodic layer. It is therefore desirable to keep the loss due 
to the resistance of the electrolyte as low as possible. 


ANODIC LAYERS 417 
happens so long as the potential of the ‘‘formed”’ 
electrode with respect to the electrolyte does not 
exceed the formation voltage on the positive 
side or from 10 to 20 percent (depending on the 
specific manner in which the layer has been 
formed and treated) of the formation voltage on 
the negative side. If the potential exceeds either 
threshold there is a considerable passage of 
current. If the positive threshold is exceeded the 
layer may form to a higher voltage or may break 
down completely, while if the negative threshold 
is exceeded the layer will become highly con- 
ducting. Therefore, if an electrolytic condenser 
is made up of a formed anode and a nonpolar- 
izable cathode and is operated on an alternating 
potential whose peak value is near that of the 
formation voltage, it will act as a rectifier. To 
prevent this, two formed anodes have to be 
used as shown in Fig. 1(A). This arrangement 
automatically biases the electrolyte at a negative 
value with respect to the average potential of 
the anodes and constitutes an a.c. electrolytic 
condenser whose capacitance is equal to the series 
capacitance of the two anodes. The potentials 
which exist in this case are shown in the oscillo- 
gram tracing of Fig. 1(A). From the picture it 
can be seen that just enough self-rectification has 
occurred so that the potential of the anode with 
respect to the electrolyte just about reaches its 
negative threshold value. This is the equilibrium 
condition, because if the threshold were exceeded, 
more rectification would take place. 

A superior operating arrangement for many 
purposes, especially those of measurement, is 
shown in Fig. 1(B). This was first proposed by 
I. Moscicki.‘ In this arrangement a potential is 
imparted to the electrolyte by separate rectifiers 
and not by the self-rectifying action of the 
anodes. Assuming that the rectifiers (4) and (5) 
rectify at zero voltage, and assuming that the 
anodes (1) and (2) have identical capacity and 
power factor, the potential of the electrolyte with 
respect to the anodes will contain a d.c. com- 
ponent of half the peak value of the alternating 
voltage impressed across the anodes. This is the 
minimum value of the d.c. component. If either 
capacity or power factor or both are different 
for the two anodes, the direct component of the 


4U.S. Patent No. 926,128. 
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(C)- ADJUSTABLE D-C ELECTROLYTE POTENTIAL 


Fic. 1. Tracings of oscillograms of electrode potentials 
as observed in a circuit (110 v, 60 cycles) with various 
methods of operation. Oscillograms: (a) é:2, potential of 
anode 2 with respect to anode 1; (b) e:3 equals potential of 
electrolyte with respect to anode 1. This is made up of 
}é:2 superimposed on the d.c. component of the electrolyte 
potential; (c) potential of anode 1; (d) d.c. component of 
electrolyte potential with respect to either anode. Circuit 
diagram: 1 and 2, formed anodes; 3, cathode; 4 and 5, 
vacuum tube rectifiers; S, sweep circuit; BT, Braun tube. 


electrolyte potential will be higher than the half 
peak value, and will increase to the full peak 
value in the extreme case when the system is 
completely unbalanced.® In the present investiga- 
tion, well-balanced (identical) anodes have been 
used in the positions (1) and (2) of Fig. 1(A), 
1(B), and 1(C). In certain special cases where 
interaction between the two formed anodes is 
suspected, critical tests have been made in which 
one anode is replaced (in the electrical circuit) by 
a mica condenser and a series resistance. 
Moscicki also suggested the introduction of an 
additional e.m.f. into his circuit which allows 


‘It is interesting to note that the minimum setting of 
the electrolyte bias in the circuit of Fig. 1(B) can be used 
for accurate capacity and power factor measurements of 
electrolytic condensers. Compare J. E. Lilienfeld, U. S. 
Patent No. 1,931,460. An analytical treatment of this 
“rectifier bridge’ is outside of the scope of the present 
paper, since the rectifier bridge was not used in any of the 
measurements reported in it. 


adjustment of the electrolyte potential to any 
desired level (Fig. 1(C)). Such an arrangement 
was used in part of the present investigation, in 
order to observe the effect of the bias of the 
electrolyte on the properties of anodes. 


TEstT SPECIMENS 


1. Superviscous specimens 


Each test specimen representing one complete 
condenser consisted of two identical separate 
cells, made in accordance with the ideas de- 
scribed in U. S. Patent No. 1,986,779. The anode 
of each cell was made of a piece of aluminum foil 
1? in. X12 in. and 10 mils thick and provided 
with a tab. The anode was formed with a direct 
current to 200 volts in all the experiments here 
described, and was inserted between two slightly 
larger copper plates with a thin coating of 
superviscous electrolyte in between. The two 
copper plates of each condenser were inter- 
connected to serve as its cathode. Two such 
identical condensers were operated in a circuit 
as shown in Fig. 1(A) in case no definite potential 
was imparted to the electrolyte, and otherwise in 
circuits as shown in Fig. 1(B) or Fig. 1(C). 


2. Aqueous specimens 


These consisted of two aluminum anodes 
1§ in.X1} in. and 10 mils thick and a copper 
cathode of similar shape in between. All three 
electrodes were mounted on a Bakelite strip that 
was supported on the edge of a container and 
were immersed in an aqueous, decidedly alkaline, 
electrolyte. The direct formation voltage of these 
samples was 155 volts. 


MEASUREMENTS 


1. General method 


We were concerned not only with capacity and 
power factor measurements, but also with a 
study of the current harmonics introduced by a 
condenser. In most of our work these current 
harmonics were so small as to be indiscernible in 
an oscillogram of the total current wave of the 
condenser; but, as we shall show later, im- 
portant information could, nevertheless, be 
learned from them. 

The arrangement sketched in Fig. 2 allowed us 
to obtain all the data presented in this paper 
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B- FOR LOW VOLTAGE EARPHONE MEASUREMENTS 


Fic. 2. Schematic diagram of measuring circuits. F, low © 


pass 60-cycle filter; O, beat frequency audio oscillator; 
PA, 250-watt, audiofrequency power amplifier; NA, null 
arm amplifier; BT, Braun tube; W, Wagner ground; C,, 
substitution resistance. 


within limits of precision adequate for our 
requirements. The equipment consisted chiefly 
of a Wien bridge whose null arm voltage was 
amplified and observed on a General Electric 
Braun tube equipped with a linear sweep circuit, 


and of oscillators covering the range from 60 to 
6000 cycles. Any harmonics present in the supply 
voltage were removed by filters. Consequently, 
all the harmonic voltages which were observed 
across the null arm of the bridge were due to 
nonlinear impedances.* The null arm amplifier 
had three resistance-capacity-coupled stages, 
having a voltage amplification adjustable from 0 
to 450. As its amplification was uniform through- 
out the range of frequencies used, the magnitudes 
of small unknown input voltages were determined 
by substituting measured 60-cycle voltages ad- 
justed to give the same deflection as the unknown 
voltage. The phase shift of the null-arm amplifier 
was constant and equal to 180 degrees from 150 
to 8000 cycles. This was calibrated by means of 
Lissajous figures on the Braun tube. Since only 
the null arm voltage is amplified in each case, 
there is no amplifier phase shift correction to be 
applied to 60-cycle voltages appearing on oscillo- 
grams. 

A great deal of the investigation was done by 
bridging over the capacity to be measured in 


®A. Gemant, “Oszillographie der Stréme in Isolier- 
stoffen,”” Archiv f. Electrotechnik 23, 683. 
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Fig. 2 with a rectifier arrangement as shown in 
Figs. 1(B) and 1(C), in order to control the 
potential of the electrolyte. For obvious reasons, 
in the range of the capacities measured (>0.1 
mf), the rectifiers should have had a negligible 
effect on the readings. This was confirmed by 
experiment. Critical tests made with good paper 
condensers, with and without the rectifier bridge, 
showed no difference in the readings to the first 
two digits even when the measured capacity 
was as small as 0.02 mf. No measurable effect 
was produced on paper condensers by varying 
the tap of the biasing battery between 0 and 
149 v d.c. Variations of the r.m.s. a.c. voltage 
between 0 and 110 v produced no changes in 
the case of paper condensers, although distinct 
variations were obtained with electrolytic con- 
densers. 


2. Capacity and power factor measurements 


To balance the bridge, when using the Braun 
tube, the circuit elements were varied until the 
fundamental disappeared on the oscillograph. 
Then, leaving the bridge settings unchanged, a 
decade mica condenser and a decade resistance 
were substituted for the test specimen and were 
adjusted until balance was again reached. The 
readings on the decade boxes then gave the 
effective capacitance and resistance of the test 
specimen. 

When using the headphones, the capacity and 
power factor of the test specimen were de- 
termined from the bridge settings. 


3. Harmonics 


When the bridge was balanced for a sinusoidal 
voltage, harmonics—chiefly the third—were in- 
troduced by the electrolytic condenser and 
appeared across the null arm.’ This null arm 
voltage, which was produced by the harmonic 
current through the condenser, amplified and 
viewed on the screen of the Braun tube, was 
photographed together with the voltage across 
the condenser and the total current through it, 
which was proportional to the drop across R, 
(Fig. 2(A)). Since the sweep frequency was 
synchronized with that of the voltage impressed 


7 No such harmonics appear with the amplification used 
if the electrolytic condenser is replaced by a good mica 
condenser. 
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Fic. 3. Phase relation of current, voltage and third 
harmonic in a superviscous condenser. A, tracing of oscillo- 
grams as taken on the circuit shown in Fig. 2. B, oscillo- 
grams A-redrawn with waves in correct phase relation. 
1, voltage applied across condenser; 2, total current 
through condenser; 3, third harmonic null arm voltage. 


on the bridge, photographs successively taken on 
the same film of voltage waves across the various 
circuit branches were shown in the same phase 
relation as if they were all taken simultaneously. 
On these photos, it had to be taken into account 
that the phase of the current through the con- 
denser was shifted by 180°. 

The impedance of the supply circuit was small 
in comparison with that of the R;C.R. branch. 
Therefore ‘the third harmonic null-arm voltage 
was in phase with the voltage drop occasioned by 
the harmonic current through R;, which was in 
phase with the harmonic current through the 
test specimen. The shifting of the switch, Sw, 
from position 1 to position 2 reverses the ap- 
parent phase of this harmonic voltage by 180°; 
but, since the null-arm amplifier has also shifted 
it by 180°, the picture of the third harmonic 
null-arm voltage shown in the photograph was 
in the same phase relation to the applied voltage 
across the condenser as the third harmonic 
current through the condenser. Thus oscillogram 
tracings such as shown in Fig. 3, reproduced 
within the limits of precision of the Braun tube, 
the correct phase relation between the applied 
voltage, the total condenser current, and the 
third harmonic condenser current, except that 
the total condenser current had to be shifted 
by 180°. 


EXPERIMENTAL RESULTS 


In Figs. 4 and 8 is shown the way in which the 
capacity and power factor of superviscous speci- 
mens vary with frequency. These variations are 
much larger than those observed in good mica 
condensers, but they are of the same order of 
magnitude and of the same general type as 
those observed in condensers with composite. 
dielectrics, such as impregnated paper. For the 
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Fic. 4. Power factor of a superviscous condenser vs. 
frequency for a constant d.c. electrolyte potential of 80 v 
and variable a.c. voltages from 0.05 to 70 v r.m.s. (at 
25°C). 80 v d.c. electrolyte potential; A, 70 v r.m.s.; B, 
50 v r.m.s.; C, 30 v r.m.s.; D, 10 v r.m.s.; E, 0.05 v r.m.s. 


convenience of practical calculations, it may be 
mentioned that variations of capacity with fre- 
quency over the range in question are closely 
approximated by the formula 


C=Co+Cie-*/ 


as demonstrated in Fig. 9. 

In the case of Fig. 4, it is more than likely 
that the part of the power loss which varies 
with frequency above 200 cycles is just that due 
to the resistance of the electrolyte. This is indi- . 
cated by the parallelity of the curves, and the 
fact that they are practically straight lines in 
this region. The smallness of this variable part 
of the power factor indicates that the total re- 
sistance of the electrolyte is very low. 

Figs. 5 and 10 show the variation of the 
capacity and power factor of superviscous con- 
densers with the d.c. component of the electro- 
lyte potential. Frequency and r.m.s. values of the 
a.c. voltage are parameters in these figures. 
The results show that the power factor varies 
but little with the biasing voltage of the electro- 
lyte.* On the other hand, the capacity does 
vary by a few percent in the permissible range of 
the electrolyte potential. All these curves appear 


8’ The biasing voltage can under no circumstances be 
allowed to exceed the formation voltage of 200 volts, 
because if it did, either the condenser would break down 
or additional layer would form and we would no longer be 
dealing with a 200-volt layer. In fact, the upper limit of the 
biasing voltage is little more than 200 volts minus one-half 
the peak of the a.c. voltage, because the layer becomes 
unstable when the formation voltage is substantially 
exceeded even in pulses. 
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Fic. 5. Power factor of a superviscous condenser vs. d.c. 
potential of electrolyte with respect to anodes. Frequency 
and a.c. voltage are parameters (at 25°C). A, 2000 cycles, 
70 v r.m.s.; B, 2000 cycles, 50 v r.m.s.; C, 2000 cycles, 
30 v r.m.s.; D, 2000 cycles, 10 v r.m.s.; E, 2000 cycles, 
0.05 v r.m.s.; W, 60 cycles, 70 v r.m.s.; X, 60 cycles, 50 v 
r.m.s.; Y, 60 cycles, 30 v r.m.s.; Z, 60 cycles, 10 v r.m.s. 
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Fic. 6. Power factor of a superviscous condenser vs. a.c. 
voltage for a d.c. electrolyte potential of 80 v and frequen- 
cies from 60 to 2000 cycles (at 25°C). 80 v d.c. electrolyte 
potential; A, 2000 cycles; B, 1000 cycles; C, 500 cycles; 
D, 240 cycles; E, 60 cycles. 


to have a characteristic minimum between 10 
and 30 percent of the formation voltage. To the 
right of these minima, the variations are approxi- 
mately controlled by a square law relation. 

Figs. 6 and 11 show how the capacity and 
power factor of a layer vary with the r.m.s. 
_ a.c. voltage for a constant electrolyte potential. 
Frequency is the parameter in these figures. 
A rising characteristic with voltage is shown in 
each case. Figs. 7 and 12 are similar to Figs. 6 
and 11, except that in obtaining the former the 
electrolyte has been allowed to establish its own 
bias. These curves are important because they 
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Fic. 7. Power factor of a superviscous condenser vs. a.c. 
voltage for frequencies from 60 to 2000 cycles with self- 
established electrolyte potential (at 25°C). A, 2000 cycles; 
B, 1000 cycles; C, 500 cycles; D, 240 cycles; E, 60 cycles. 
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Fic. 8. Capacity of a superviscous condenser vs. fre- 
quency at fixed d.c. electrolyte potentials with a.c. voltage 
as parameter (at 25°C). A, 110 v r.m.s., 77.5 v d.c. electro- 
lyte potential (3 a.c. peak voltage); B, 70 v r.m.s., 80 v 
d.c. electrolyte potential; C, 60 v r.m.s., 80 v d.c. electro- 
lyte potential; D, 40 v r.m.s., 80 v d.c. electrolyte potential; 
E, 20 v r.m.s., 80 v d.c. electrolyte potential; F, 0.05 v 
r.m.s., 80 v d.c. electrolyte potential. 


picture the conditions under which practical 
condensers are usually operated. 

It should be mentioned that the variations of 
the capacity and power factor with frequency 
are of a nonlagging type, but the variations of 
these with the electrolyte potential and the 
r.m.s. a.c. voltage are of a lagging character. 

In Figs. 13-15 and Fig. 3, a study is made of 
the harmonics generated in electrolytic con- 
densers. In Fig. 13 are shown oscillogram tracings 
at different frequencies, each showing the applied 
potential, the total condenser current and the 
amplified null-arm voltage. The latter gives a 
picture of the harmonics generated by the con- 
denser. At 60 cycles there is practically a pure 
third harmonic. The third harmonic remains the 
preponderant component at all the frequencies 
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Fic. 9. Plot showing that the capacity per unit area, C, 
varies with the frequency, f, according to the relation, 
C=(Co+C,%. A, capacity per unit area vs. frequency; 
B, log (C—Co) X10* vs. frequency. C=capacity per unit 
area as a function of frequency. Co=capacity_at infinite 
frequency. 
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Fic. 10. Capacity per unit area of a superviscous con- 
denser vs. d.c. electrolyte potential, with frequency and 
a.c. voltage as parameters (at 25°C). A, 70 v r.m.s.; B, 
50 v r.m.s.; C, 30 v r.m.s.; D, 10 v r.m.s.; E, 0.05 v r.m.s. 
Circles with vertical line, 60 cycles; circles with horizontal 
line, 500 cycles; circles with lower half solid, 2000 cycles. 


photographed, but other harmonics become more 
and more in evidence as the frequency is in- 
creased. By taking into account the amplification 
of the amplifier, measurements show that the 
amplitude of the harmonic current is roughly 
the same at all the frequencies in question. The 
curves of the total condenser current show that 
the harmonic current is very small in absolute 
magnitude, in fact of the same order as that ob- 
served in good impregnated paper condensers. 
In Fig. 14 is shown that the amplitude of the 
third harmonic current is proportional to the 
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Fic. 11. Capacity of a superviscous condenser vs. a.c. 
voltage, for a d.c. electrolyte potential of 80 v with fre- 
quency as parameter (at 25°C). 


cube of the applied voltage at 60 cycles. A similar 
relation holds at 240 cycles, and seems to hold 
even at higher frequencies, to the limited extent 
to which measurements can be made. 

The percentage of harmonics varies with the 
electrolyte in which the layer is formed and 
operated. By using aqueous specimens with a wide 
range of OH~ concentrations, Fig. 15 is obtained. 
It shows that the percentage of harmonics is 
roughly proportional to the minimum formation 
current of the electrolyte.® 

In Fig. 3 (A and B) are shown the phase rela- 
tions between the voltage, total condenser 
current, and harmonic current in a superviscous 
specimen at 60 cycles. To a close approximation, 
the corrected curves of Fig. 3(B) indicate that 


e=E sin wt, 
i=I, cos wt—Ts sin 3wt, 


where e and 7 are the total voltage and current, 


® In a report to be published in Trans. Electrochem. Soc. 
October (1935), it is shown that for any ‘‘forming’”’ electro- 
lyte in which an aluminum anode is immersed, there is a 
critical current density at the Al anode above which for- 
mation of the anode takes place. The value of this critical 
current density, which is called the minimum formation 
current, is independent of the voltage to which the anode 
has already been formed. The minimum formation current 
is approximately proportional to the OH~ concentration of 
the electrolyte. 
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Fic. 12. Capacity vs. a.c. voltage for frequencies with 
self-maintained electrolyte potential (at 25°C). A, 60 
cycles; B, 240 cycles; C, 500 cycles; D, 1000 cycles; E, 
2000 cycles. 
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Fic. 13. Tracings of photographs showing the properties 
of the harmonic current through the condensers from 60 to 
1000 cycles. 1, amplified null arm voltage (representing 
harmonic current); 2, voltage of electrolyte with respect 
to anodes; 3, total current through the condenser. 


respectively. Similar relations hold at higher 
frequencies, except for the increasing evidence 
of other harmonics. 


THEORETICAL CONSIDERATIONS 


Nonlagging effects 


The frequency characteristics of capacity and 
power factor indicate that the dielectric of an 
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Fic. 14. Variation of third harmonic current and its 


cube root vs. applied voltage at 60 cycles for a superviscous 
condenser. 
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the minimum formation current of the electrolyte (110 v 
r.m.s. at 60 cycles across the condenser). 
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anodic layer is similar in properties to other 
composite dielectrics. Its characteristics fit in 
quite well with Maxwell’s absorption theory as 
extended by K. W. Wagner." The composite 
nature of the dielectric of an anodic layer is 
probably due to the occlusion of small amounts 


10 Compare F. W. Grover, Bull. Bur. Standards 7, 495 
(1911). 

J. B. Whitehead, Lectures on Dielectric Theory and 
Insulation (McGraw Hill, New York, 1927)—Lectures 6 
and 7. 
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of aluminum hydroxide in the matrix of alumi- 
num oxide. 

In the section on harmonics, it will be shown 
that there is a nonlagging variation with voltage 
of the a.c. conductivity. 


Lagging effects 


The characteristics of anodic layers also 
present a class of phenomena peculiar to these 
layers and caused by their dynamic nature. 
These phenomena are lagging effects and indicate 
changes in the organization of the layer. They 
include the variation with a.c. and with d.c. 
voltage of the capacity and power factor. A de- 
tailed interpretation of these effects would be 
very difficult, pending a better knowledge of the 
ultimate factors which govern the capacity and 
power factor of anodic layers. 


Harmonics 


A diagram of a network equivalent to an 
anodic layer in a.c. operation is shown in Fig. 16, 
where C represents the capacity; R, the power 
factor; and G, the a.c. leakage conductance. 
Now, the presence of a preponderantly third har- 
monic in the harmonic current must be accounted 
for by a quadratic variation with voltage in one 
of these elements. A comparison of the experi- 
mental results with the results to be obtained by 
a variation of each of these factors separately, 
indicates that the quadratic variation occurs in 
the a.c. leakage conductance. 


Let e=E sin wt (1) 


be the variable part of the e.m.f. impressed on 
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Fic. 16. Equivalent circuit of an electrolytic condenser 
when operated on an a.c. circuit. C, equivalent capacity 
of electrolytic condenser; R, equivalent series resistance 
representing power factor; G, equivalent leakage con- 
ductance to alternating current; Go, limiting value of G 
for small a.c. voltages; K, proportionality constant; 7, 
total current through condenser; ic, current through 


equivalent capacity; tg, current through equivalent a.c. 
leakage conductance. 


e=E sinwt G=Gi+Ke) 
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the circuit of Fig. 16. Assume 
G=G,(1 + ke’), (2) 


where Gp may be a function of the biasing 
voltage.” Then 


EwC 
=———— cos wi+ 
1+ R*w*C? 


ER«*C? 


in wt 


es Si 
1+ R*w*C? 


3k 
+G E(1+ 4k’) sin - Gr sin 3wt. (3) 


Then, realizing that the second and third terms 
of the expression on the right are small in com- 
parison with the first, we see that the above gives 
the correct phase relation between the harmonic 
current and total current. (Compare Fig. 3.) 
It likewise shows that the amplitude of the har- 
monic current should vary as the cube of the 
impressed a.c. voltage (compare Fig. 14) and 
accounts for the fact that the amplitude of the 
third harmonic current varies little, if at all, 
with frequency (Fig. 13). 

A corresponding analysis with a wonlvatle 
variation of the capacity, instead of the a.c. 
leakage conductance, gives the wrong phase for 
the harmonic current; and also requires that it 
increase with frequency, which is likewise not 
observed. The analysis with a quadratic variation 
of the power factor, while giving the correct 
phase relations, requires that the amplitude of 
the harmonic current rise wjth frequency, in 
contradiction to the experimental results. We are 
thus led to the conclusion that the third harmonic 
current is due to a quadratic variation in the a.c. 
leakage conductance. This conclusion receives 
additional support from Fig. 15, where it is 
shown that the percentage of harmonics is 
directly proportional to the minimum formation 


12 This expression, G=G»(1+ke*), means that the total 
conductivity is of the form 


& = Goler){ 1+k(e’ —e1)*} +f(e1) (A) 


where e’ is the total potential drop across the layer and 
é, is the biasing voltage of the electrolyte. Some experi- 
ments which we have performed indicate that & is directly 
proportional to the biasing voltage. This variation of k 
with the biasing voltage is, however, a lagging a. - 
it is not possible to express it explicitly i in Eq. (A). T 
significance of Eq. (A) is that with each new biasing ae 
age, the anodic layer passes into a new equilibrium 
condition for which a current-voltage relation of the form, 
Eq. (2), holds. 
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current.” The agreement of the calculations with 
the results observed in Figs. 3 and 13 indicates 
that the time lag, if any, between the change of 
voltage and the change of conductance is less 
than 10-* second. 

It is not to be expected that the relation 
between leakage current and voltage derived 
from the presence of the third harmonic can be 
confirmed by reading the d.c. leak at consecutive 
voltages with a d.c. instrument. In the latter 
procedure, lagging changes of leak are super- 
imposed upon the nonlagging changes and each 
point is unavoidably taken with the layer in a 
changed condition. The result of such a slow 
observation will vary according to the procedure 
of the experiment, i.e., how fast it is done, 
whether voltages are increased or decreased, and 
other arbitrary factors. Therefore, d.c. current 
vs. voltage curves have very limited significance 
for anodic layers. On the other hand, the results 
of our experiments on harmonics show that (in 
carefully made-up superviscous condensers and 

13 Although it is not possible to get a direct measurement 
of the leakage current of an anodic layer in a.c. operation, 
it is known that the d.c. leak is about proportional to the 


minimum formation current, and it is therefore reasonable 


to suppose that in a general way a similar rule holds for 
the a.c. leak. 
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in such aqueous condensers as are not too 
alkaline and do not contain detrimental ions) for 
instantaneous variations of current with voltage, 
the relation 1=G)(1+ke?) holds. 

The absence of harmonics of such magnitude 
as would correspond to the variations of capacity 
with voltage shown in Figs. 10 and 11, isadditional 
evidence that these variations are lagging effects. 
It should be mentioned that a number of in- 
vestigators have looked for capacity changes 
that could be correlated with the compressibility 
of the layer under the mechanical pressure on it 
caused by the electrostatic field. This pressure 
reaches in the neighborhood of 1000 kg/cm’. 
Now, capacity changes due to pressure would be 
of a ‘“‘nonlagging”’ character, and the oscillo- 
graphic investigation of harmonics used here 
should bring them to light. Therefore, the 
absence of evidence of such changes of capacity 
indicates that the compressibility of the layer is 
very small. From this we may conclude that the 
layer cannot be of a gaseous nature, as has been 
assumed by some earlier investigators. 

The authors wish to mention that Mr. Alfred 
Herzenberg was connected with this work during 
its initial period. 
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